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CHAPTER 1 
One-Dimensional Quantum Spin Liquids 
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Universite de Cergy-Pontoise, 5 mail Gay-Lussac, Neuville sur Oise, 
95301 Cergy-Pontoise Cedex, France 



This chapter is intended as a brief overview of some of the quantum spin 
liquid phases with unbroken SU(2) spin symmetry available in one di- 
mension. The main characteristics of these phases are discussed by means 
of the bosonization approach. A special emphasis is laid on the interplay 
between frustration and quantum fluctuations in one dimension. 



1. Introduction 

A central issue in the study of strongly correlated systems is the classifica- 
tion of all possible Mott insulating phases at zero temperature. The general 
strategy for describing the possible phases associated with the spin degrees 
of freedom (the so-called "quantum magnetism" ) of Mott insulators con- 
sists of analysing localized spin models such as the antiferromagnetic (AF) 
Heisenberg model: 

U^jY. S,-S, + .., (1) 

<ij> 

where Si are quantum spin-S" operators on sites i of a lattice, and the sum 
< i,j > is over the nearest neighbor sites with an antiferromagnetic ex- 
change ( J > 0); the ellipsis represents additional terms like second- neighbor 
competing exchange interaction or ring exchange that can be eventually 
added to the Heisenberg Hamiltonian (Q. Generically, the ground state of 
the model without these additional terms displays long-range Neel or- 
dering (< Si >^ 0) which breaks spontaneously the SU(2) symmetry of the 
lattice Hamiltonian (0). The low-energy excitations are gapless spin- waves 
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(i.e. magnons) as expected when a continuous symmetry is spontaneously 
broken. 

A central focus of quantum magnetism over the years, starting from the 
proposal made by Anderson^, has been the search for a spin liquid behavior 
i.e. a phase with no magnetic long-range Neel order. Spin liquids phases are 
expected to be stabilized in low dimensions or in presence of frustration, i.e. 
in situations where quantum fluctuations can strongly suppress magnetism. 
It gives rise to rich physics with exotic low-energy excitations which require 
in most cases the use of non-perturbative techniques to fully determine its 
properties. In the two-dimensional case, different spin liquid phases have 
been found and display bond ordering or topological ordering (for a review 
see for instance Refs. "^Ml and references therein). An interesting attemp to 
develop a systematic classification of two-dimensional spin liquids has also 
recently been explored by Wen El. 

In one dimension, the quantum analog of the Mermin- Wagner El theo- 
rem shows that quantum fluctuations always disorder the ground state in 
systems with a continuous symmetry. The generic situation in one dimen- 
sion is thus to have a spin liquid phase. The one dimensional case is also 
extremely favorable since several non-perturbative techniques are available 
to fully characterize the physical properties of different spin liquids. T hese 

deintegr; 

3i]m, 

density matrix renormalization group (DMRG) approach ^5 The determi- 
nation of the properties of spin liquids is not a purely academic problem 
since many physical realizations of spin chains have been synthes ized ov er 
the years. The reader may consult for instance the recent reviews ^^1 on 
experiments in quasi-one dimensional spin systems. 

The main questions in the study of one-dimensional spin liquids are 
in order: What are the different spin liquid phases available at zero- 
temperature? What are the different physical properties of these spin liq- 
uids? What is the nature of the quantum phase transition 13 between two 
different spin liquids? Last but not least, the problem of dimensional cross- 
over, i.e. the fate of the spin liquid as the two dimensional case is reached 
is also a central topic. In particular, is it possible to stabilize a spin liquid 
phase in d = 2 starting from the d — 1 case? 

In this work, we shall mainly be concerned with the two first questions 
for spin liquids with unbroken SU(2) spin symmetry. Using the bosonization 
approach, we shall discuss some of the physical properties of spin liquids 
that have been found in spin chains, spin ladders with or without frustra- 



powerful techniques includ e integ rability, conformal field theory (CFT)I^EI 
the bosonization approacl0^^^^ and numerical calculations such as the 
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One- Dimensional Quantum Spin Liquids 3 

tion. An important question related to this problem is how can one clas- 
sify all these ID spin liquid phases at T = 0? A first distinction concerns 
the existence or absence of a spectral gap in the model. In a gapless sys- 
tem (critical spin liquid), one has quasi-long range Neel order. The leading 
asymptotics of spin-spin correlation functions display a power law behavior 
with an exponent characterized by a well defined underlying CFT ■ . In 
contrast, in a spin liquid phase with a spectral gap, spin-spin correlation 
functions have an exponential decay with the distance due to the existence 
of a finite correlation length ^: {Si ■ Sj) ^ exp(— |i — j\/£,). At this point, 
it is worth stressing that two gapful spin liquid phases sharing the same 
thermodynamic properties (typically a thermal activation law at low T) 
may display very different behaviors in other physical quantities of interest 
such as, for instance, dynamical or optical properties. It is thus necessary to 
make a scrutiny analysis of the properties of the phase before elaborating 
a complete classification. Let us first consider the ground state of a gapful 
spin liquid phase. A (discrete) symmetry might be spontaneously broken 
in the ground state resulting on a ground-state degeneracy. The spin liq- 
uid phase may display also an hidden topological order ^1 which manifests 
itself in a ground-state degeneracy which depends on the nature of bound- 
ary conditions (BC) used i.e. periodic or open BC. The resulting chain-end 
(edge) excitations in the open BC case can be observed in the NMR pro- 
file of a spin chain compound doped with non-magnetic impurities like Zn 
or Mg. Another important distinction between two spin liquids relies on 
the quantum number carried by elementary excitations. The fractionalized 
or integer nature of this quantum number has crucial consequences in the 
dynamical structure factor of the system which directly probes the elemen- 
tary nature of the spin-flip. A sharp spectral peak, characteristic of a well 
defined = 1 mode, or a broad (incoherent background) feature in the 
dynamical structure factor can be seen in inelastic neutron scattering ex- 
periments. Finally, an interesting question is the existence or not of bound 
states in the energy spectrum below the two-particle continuum. Spin liquid 
phases may then be distinguished at this level which is also an important 
fact from the experimental point of view since observation of magnetic sin- 
glet bound states can be realized through light scattering experiments 
This classification is certainly far from being complete but it enables us to 
distinguish between several spin liquid phases and investigate the nature of 
the quantum phase transition between them. 

This chapter is organized as follows: We present, in Section 2, the dif- 
ferent spin liquid phases that occur in unfrustrated spin chains and spin 
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ladders. The main effects of frustration in one-dimensional spin liquids are 
described in Section 3. In particular, we shall observe that frustration plays 
its trick by allowing deconfined spinous (carrying fractional S* = 1/2 quan- 
tum number) as elementary excitations and it provides a non-trivial source 
of incommensurability. Finally, our concluding remarks are given in Section 
4. 



2. Unfrustrated spin chains 

The paradigmatic model to investigate the properties of ID quantum mag- 
nets in absence of frustration is the AF Heisenberg spin chain given by the 
Hamiltonian: 

n = JY,S^■S^+u (2) 

i 

where Si is a spin-S* operator at the ith site of the chain and the exchange 
interaction is antiferromagnetic: J > 0. In the classical limit, the model 
displays Neel long-range order: Si ^ S'(— l)'r?, n being an unit vector with 
an arbitrary fixed orientation in spin space. This solution breaks the SU(2) 
symmetry of the model Q down to U(l). Corresponding to this symme- 
try breaking scheme, there are two Goldstone modes which propagate with 
the same velocity (Lorentz invariance in the low-energy limit). In the lan- 
guage of CFT 13 it corresponds to a field theory with central charge c = 2 
{c — I being the central charge of a free massless boson field). Physically, 
these Goldstone modes are nothing but a doublet of gapless spin waves 
modes associated with slow modulations in the orientation of the vector ft 
which represents the order parameter of the Neel magnetic structure. In 
the quantum case, this Neel solution in one dimension is destabilized by 
strong quantum fluctuations no matter how large S is: a spin liquid phase 
is formed. 



2.1. Spin- 1/2 Heisenberg chain 

In the ultra quantum limit i.e. S = 1/2, the nature of the spin liquid phase 
can be fully determined since the model is exactly solvable for S* = 1/2 
by means of the Bethe ansatz approach. In particular, starting with Bethe's 
seminal work^l, a host of exact results have been ob tain ed over the years 
for ground state properties'^ magnetic susceptibility^ ^ th ermodynamics 
123 excitation spectrum 1^22^ and correlation functions 

The model displays quantum criticality properties which belo ng t o the 
Wess-Zumino-Novikov-Witten (WZNW) su(2)i universality class l^^. The 
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One- Dimensional Quantum Spin Liquids 5 

low-energy limit of the spin- 1/2 AF Heisenberg chain Q is described by 
the su(2)i WZNW CFT with central charg e c = 1 perturbed by marginally 
irrelevant current-current interactional^. The resulting Hamiltonian den- 
sity reads as follows: 

2t:v 



He 



^ {Jl + 4) + A/fi • Jl, (3) 



where Jr and Jl are respectivel y th e right and left su(2)i currents which 
generate the su(2)i WZNW CFtES these currents satisfy the su(2)i Kac- 
Moody commutation relations: 



[JIl (^) , Jr.l (y)] = {x~y) + ze'^^V]^,^ {x) 5 (x - y) . (4) 

In Eq. , w is the spin velocity and A < so that the last contribution is a 
marginally irrelevant term that renormalizes to zero in the far infrared (IR) 
limit. This pert urb ation accounts for logarithmic corrections in the spin- 
spin correlation 1^ which is exactly known in the long-distance limit: 

The most striking feature of this spin liquid phase stems from the nature 
of its elementary excitations. They have been elucidated by Faddeev and 

on 

Takhtajan within the Bethe ansatz approach and consist of fractional 
= 1/2 massless excitations called spinous. The lowest excitations are 
fourthfold degenerate and correspond to a triplet (S* = 1) and a singlet 
{S = 0). The resulting energy spectrum is a continuum in (k,uj) space 
bet ween a lower boundary (the des Cloizeaux-Pearson dispersion relation 
l^^l) Wdcp = 7rJ|sinfc|/2 (— tt < k < it) and an upper boundary uju = 
7rJ| sin(fc/2)| (see Fig. The central point of the analysis is that this 
continuum can be interpreted as being made up of two spin-1/2 excitations 
(spinous) with the dispersion: 

ttJ 
"2 

with < fc < TT. A spinon has thus a wave-vector restricted to only half 
of the Brillouin zone. A triplet (or singlet) excitation with momentum k is 
then described by two spinous with momenta ki and ^2 (0 < fci < ^2 < tt) 
such that: uj{k) — Wspinon(^i) + '^spinon(fc2) with fc = fci -|- ^2 (respectively 
k = ki + k2 — 2tt) if < fc < tt (respectively — vr < fc < 0). It corresponds 
to a two-parameter continuum uj{k) = 7rJsin(fc/2) cos(fc/2 — qi) with < 



spinon — ~;^sinfc, (6) 
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gi < k/2 for < fc < tt and 7r + fc<gi <7r + fc/2 for — tt < < which 
identifies with the continuum of the spin-1/2 AF Heisenberg chain. 

(OA 




Fig. 1. Two-spinon continuum of the spin-1/2 AF Heisenberg chain. 

A magnon excitation carrying S = 1 quantum number is thus no longer 
an elementary excitation in this model but is fractionalized into two spinous. 
One important experimental signature of these spinous is the absence of any 
sharp peak in the dynamical susceptibility structure factor since a single 
spin-flip generates a triplet excitation which is not elementary here but 
made up of two spinous. The spectral density is then a convolution of the 
spectral densities of the individual spinous. These spinous have fascinating 
properties in particu lar they obey semion statistics intermediate between 
bosons and fermions as it can be seen by the study of the elementary 
excitations of the Haldane-Shastry model *^ which belongs to the same 
universality class as the spin-1/2 AF Heisenberg chain. An heuristic way 
to describe the spinous is to consider an anisotropic XXZ version of the 
Heisenberg model (jJJ: 

Wxxz ^ Jy-'^ '5'f <5'f+i + X! {^t^'i+i + ^t+i^i ) ■ C^) 

i i 

This model is still integrable and the ground state is two-fold degenerate for 
Jz > Ji_ as the result of the spontaneous breaking of a discrete Z2 symme- 
try. In this regime, a spectral gap is opened and the elementary excitations 
are massive spinons. A simple way to understand these properties is to con- 
sider the Ising limit when Jz ^ J± where the ground state reduces to the 
two Neel ordered states: |J,|J,|J,-- and J,|J,|J,|--- The elementary exci- 
tations are massive kinks which carry AS^ — ±1/2 quantum number and 
interpolate between these two ground states: |J,|J,| x |J,|.. The compos- 
ite nature of a magnon excitation can then be readily seen by starting from 
a ground state and apply a single spin-flip on it to obtain the configuration: 
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One- Dimensional Quantum Spin Liquids 7 

TJ-imtJ-Ti' ^PPly™S tti^ exchange J± interaction of the model (0 on 
this state, this magnon excitation with AS^ = —1 is in fact made up of two 
kinks of the Neel order i.e. two spinons: ||x||J,|J,xJ,|J,-At the transi- 
tion Jz ~ J±_ (isotropic Umit) , the ground state is disordered by "condensa- 
tion" of these kinks and the low- lying excitations carrying AS"^ = ±1/2 are 
transmuted from gapful kinks to gapless spinons. Magnons are still built as 
pairs of these spinons and have a gapless energy spectrum. 

This remarkable production of two particles from a single spin-flip has 
been explored in inelastic neutron scattering measurements on quasi-lD 
compounds that are good experimental realizations of the model The 
first clear evidence of spinons comes from neutron experiments^^on KCuFa 
that found a continuum of magnetic excitations consistent with that ex- 
pected from unbound spinons pair. The spinon continuum has also been 
observed in several other realizations of the spin - 1 1 2 AF Heisenb erg c hain 
as the copper benzoate Cu(C6D5COO)2.3D2 ED, BaCuaSiaOr ^S, and 

Cu(C4H4N2)(N03)2'^. 

2.2. Haldane's conjecture 

After the elucidation of the elementary excitations of the spin- 1/2 AF 
Heisenberg chain by Faddeev and Takhtajan E21 in the early eighties, a 
central question of ID quantum magnetism concerns the stability of the 
spinons. What is the fate of these exotic excitations upon switching on 
small deviations from the model (01? Do these additional terms lead to a 
confinement of the spinons or do they still remain deconfined? 

A simple modification of the spin-1/2 AF Heisenberg chain (jSJ is not to 
add perturbations but to change the Hilbert space namely to consider the 
general spin-S* case. The resulting model is no longer integrable and one 
has to resort to approximate methods to describe its physical properties. 
This question has lead to one of the most profound result in the field of ID 
quantum magnetism: the difference between i nteg er and half-integer spins 
AF Heisenberg chains predicted by Haldane 1"^. Integer spin chains are 
incompressible spin liquids with a finite gap (the so-called Haldane gap) 
in the energy spectrum. In that case, the spinons are confined and the 
elementary excitations of the model are massive (optical) 5* = 1 magnons. 
In contrast, for half-integer spins, the model displays quantum criticality 
with similar universal properties as the 5=1/2 case that are described by 
the su(2)i WZNW universality class. The elementary excitations are still, 
in this case, gapless spinons. 
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This remarkable distinction between integer and half-integer spins, his- 
torically called Haldane's conjecture, has been obtained by means of a semi- 
classical analysis of the model with a special emphasis on the topological 
nature of the order parameter fluct ua,tions . For a review of this approach, 
the reader may consult the reviews l^'^ l -^^ l and the books l36|(37,14| 
large S limit, the Euclidean action that describes the low-energy properties 
of the spin-5 AF Heisenberg chain reads as follows 



where v = 2JS, g = 2/S, and n is the order parameter of the Neel 
collinear state. The second contribution in Eq. ||SJ) is a topological term 
since Q {ft) — J drdx n- {drfi /\dxn) / An is an integer, called the Pontryagin 
index, which measures the number of times the spin configuration fi{x, r) 
covers the surface of the unit sphere S^. In more mathematical terms, the 
configurations r), with fixed boundary conditions at infinity, are map- 
pings of the sphere onto with homotopy classes classified by an integer 
112 ('5'^) — Z , which is nothing but the Pontryagin index Q. For integer spin 
chains, the term 2ttSQ (n) in Eq. © has no effect on the path integral 
of the model and can be discarded. The effective action ^ reduces then 
to the one of the two-dimension al 0(3 ) non-linear sigma model which is 
a massive integrable field theory popi^ l xhe exact low-energy spectrum of 
this field theory consists of a massive bosonic triplet with mass m ~ e"^'^/^ 
and it exhibits no bound states. The AF Heisenberg chain with integer spin 
is thus expected to have low-lying triplet excitations with a gap that scales 
as A ~ ve~'^^ when 5* is large. 

In the case of half- integer spins, the topological term in Eq. mani- 
fests itself in the path integral through a phase factor (—1)^ which gives rise 
to quantum interference between topologically distinct paths in space-time 
of the order parameter field n{x,T). It turns out that this process protects 
the model from a dynamically generated mass gap by quantum fluctuations 
and, in contrast, it leads to a non-p erturbati ve massless flow towards an 
IR conformally invariant fixed point EH^^SSl -vvhich belongs to the su(2)i 
WZNW universality class. Since this fixed point describes also the univer- 
sal properties of the spin-1/2 AF Heisenberg chain, all half-integer spins 
chains should exhibit the same type of emerging quantum criticality with 
for instance spin-spin correlations that decay according to the power law 
behavior ij^J. 




(8) 
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2.3. Haldane spin liquid: spin-1 Heisenberg chain 

The simplest incompressible one-dimensional spin liquid phase corresponds 
to the spin-1 AF Heisenberg chain. The existence of a Haldane gap in this 
model has generated an intense activity over the years after its prediction. 
It has been confirmed numerically from exact diagonalizations on finite 
samples 133, quantum Monte-Carlo methods Ei, transfert matrix compu- 
tationsl^, and finally from DMRG calculations In particular, this 

latter technique predicts a gap A = 0.41050(2) J and a correlation length 
^ ~ 6.03(1) lattice spacings. From the experimental point of view, several 
Hal dane compounds have been synthesized over the years (see for instance 
Ref. for a recent review) . The two most studied compounds are CsNiCla 
E3and Ni(C2H8N2)2N02(C104) (NENP) ESI where spin-1 local moments 
are provided by Ni^+ ions. In particular, inelastic neutron scattering ex- 
periments on NENP confirm the existence of the Haldane gap. This 
energy gap has also b een observed i n se veral other q uas i- 1 D sp in- 1 materials 
as (CH3)4NNi(N02)3Eni, AgVP2S6EIl, and Y2BaNi05E3Slwhich is prob- 
ably the best realization of the spin-1 AF Heisenberg chain. Finally, it is 
worth noting that the recent compounds Ni(CaDi4N2)2N3(PF6) (NDMAP) 
Eland Ni(C5Hi4N2)2N3(C104) (NDMAZ) enable to investigate exper- 
imentally the high magnetic field properties of the spin-1 Heisenberg chain 

From the theoretical point of view, several methods have been intro- 
duced to shed light on the nature of the mechanism of the Haldane-gap 
phenomena and to determine the main characteristics of this incompress- 
ible spin liquid phase. These approaches ar e the non-linear sigma model 
field theory obtained in the large spin limit l^^ l ^^ l the valence bond state 
(VBS) description!^, t he Mai orana fermions methodic, and the restricted 
Hilbert space approach 

Probably, the simplest and most appealing approach to study the phys- 
ical properties of the spin-1 AF Heisenberg chain consists to add a bi- 
quadratic interaction to the model 



One of the main interest of this extended Heisenberg model is that for 
f3 = 1/3, the so-called AKLT point ^' , the model has an exactly solvable 
ground state which captures the main characteristics of the Haldane spin 
liquid phase. At this special point, the Hamiltonian jnj is equivalent to 
the sum of projection operators P2{i,i + 1) that project onto the spin-2 
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contribution for every pair of nearest-neighbor spins: 



KLT 



n 



/9=l/3 



2J^ 



P2(^,^ + l) 



(10) 



From this structure, the ground state of this model can be constructed ex- 
actly using nearest-neighbor valence bonds. Following Affleck and coworkers 
pionneer workl^, each original 5 = 1 spin is written as two spin-1/2 vari- 
ables in a triplet state. The ground state is then obtained by coupling into 
a singlet state all nearest-neighbor spin-1/2, thus forming a crystalline pat- 
tern of valence bonds (see Fig.EJ. This state is called the valence-bond-solid 




Fig. 2. VBS state in a chain with periodic boundary conditions. 

(VBS) state. For periodic BC and an even number of sites, this VBS state is 
the unique ground state and does not break the tr ansl ation symmetry of the 
Hamiltonian H10() . It can also be shown rigorously!^ that the AKLT model 
has a spectral gap whose value has been computed numerically: A ~ 0.7 J 

In addition, the authors of Ref. have also been able to show that 
the spin-spin correlation function decays exponentially with the distance 
with a finite correlation length equals to ^ = l/ln3 lattice spacings. The 
VBS state is therefore the simplest incompressible spin liquid phase with a 
unique ground state and a gap to all excitations. 

Despite the exponential decay of correlations in this spin liquid phase, 
there is a subtle form of hidden AF orde ring. This hidden order was first 
discovered by den Nijs and Rommelse using an equivalence of the spin-1 
chain t o a two-dimensional restricted solid-on-solid model and also later by 
TasakiESlby means of a geometrical approach. A simple way to exhibit this 
hidden AF ordering is to look at the VBS state written in the conventional 

representation. A typical configuration has the following structure: ..0 t 
0..0 itit i 0..0 t 0.., i.e. each t (S^ = +1 state) is followed by i (fif, = -1 
state) with an arbitrary number of states (5^ = state) between and vice 
versa. If we ignore the spins with = state then the remaining spins 
display a perfect long-range spin-1/2 Neel ordering: tititi •■• The VBS 
state has thus a perfect dilute AF Neel order. Because of the arbitrary 
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One- Dimensional Quantum Spin Liquids 11 

number of S*^ = sites inserted between Sp — ±1 sites, the long-range 
AF order is invisible in the spin-spin correlation function which has an 
exponential decay. However, this hidden order becomes manifest in the 
non-local string order parameter introduced by den Nijs and Rommelse 1^ 
and defined by 



a^i^g = - hm (Sfexp ( ^ | Sf), (11) 



with a = x,y, z. In t he VBS ground-state, this order parameter can be 
computed exactly O"^^.^^^ = 4/9 i.e. this object exhibits long-range 
order in the VBS phase. The existence of such a hidden order is essential 
to the basic mechanism of the Haldane gap: breaking this topological order 
to create an excitation costs a finite energy gap. The properti es o f this 
dilute AF order was further elucidated by Kennedy and Tasakil^. They 
were able to show, using a non-local unitary transformation, that the long- 
range order Ofj^j^g ^ and the Haldane gap are related to a spontaneous 
breaking of a hidden Z2 x Z2 symmetry. In particular, the string order 
parameter becomes the usual local ferromagnetic order parameter under 
this non-local unitary transformation. Another important consequence of 
this broken Z2 x Z2 symmetry is the existence of a quasi-degeneracy of 
four lowest energy levels for a finite-size chain. This can also be understood 
within the VSB state since for open BC the AKLT model (|10|l has exactly 
four ground states with the existence of two spin-1/2 degrees of freedom 
that are unpaired at each end of the chain (see Fig. EJ. Three of these 







Fig. 3. VBS state with open boundary conditions; unpaired bonds are left at the bound- 
aries resulting into two free spin-1/2 objects at the edges. 



states constitute a spin triplet whereas the last one is a spin singlet. 

Beside the fact that only the ground state of the Hamiltonian H10|l is 
known exactly, the AKLT approach is also very useful to describe low-lying 
excitations of the model throujgh a variational approach starting from the 
VBS state. A first approach E21 consists to create a "magnon" excitation by 
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applying a spin-flip on the VBS ground state {\VBS)): 




(12) 



N being the system size. It is very tempting to interpret this excitation 
as a standard S = 1 magnon i.e. an excitation which modifies only locally 
the AF Neel order. In fact, the trial wave function 1)12(1 does not describe 
a magnon excitation of the perfect hidden Neel order of the VBS state but 
it has instead a solitonic nature Indeed, in the representation, a 
typical configuration that appears in Sj'\VBS) has the following structure: 
..0 t 0..0 it t i 0..0 t 0.. Removing the states, this configuration 
becomes tittit ■• ^ soliton or domain wall which destroys non-locally 
the AF Neel long-range order. An alternative approach to describe the 
elementary excitations of the AKLT Hl()|l model is to start from the VBS 
state and to promote the spins of the link + 1) into a triplet state |$°) 
I ti^Kjl l ^? ] This low-lying excitation carrying momenta k describes a moving 
triplet bond, usually called crackion, and it is defined by 



In fact, the two approaches (|12|l and H13() are related and describe the same 
dispersion relation for the triplet excitation The solitonic nature of 
the excitation (|13() becomes manifest by considering open boundary con- 
ditions an d th e non-local unitary transformation introduced by Kennedy 
and TasakiEl. In that case, the trial excited wave function corresponds to 
a domain wall that interpolates between two of the four possible ground 
states 

One of the main success of this VBS approach stems from the fact that 
the properties of the VBS state are not special to the AKLT point ((l()|l but 
have a larger extent. This model turns out to be smoothly connected to 
the spin-1 Heisenberg chain in the sense that they s hare same physic al 
properties. In this respect, numerical investigations I v^lv'^l < v l Y 1| (A (41 1 4 | qf 
the phase diagram of the bilinear-biquadratic spin-1 Hamiltonian Q have 
revealed that this model with \(3\ < 1 belongs to a same phase, the Haldane 
phase, whose main characteristics are well described by the VBS state. In 
particular, the origin of the Haldane gap in the spin-1 Heisenberg chain is a 
consequence of a hidden Z2 x Z2 broken symmetry as for the AKLT model 
((l()|l . At the Heisenberg point f3 — 0, the hidden AF Neel ordering is no 




(13) 
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longer perfect as at the AKLT point. It is only locally destroyed so that 
there is still long-range order at /? = with a non-zero string ord er par arn- 
eter (HTJ but smaller than at the AKLT point: O^^Hng - 0.374 llUiZlIlSl. 
A second consequence of the spontaneous breaking of the Z2 x Z2 sym- 
metry is the two almost free S = 1/2 degrees of freedom at the ends of a 
finite chain with open BC. At the Heisenberg point, exact diagonalizations 
1221 of finite open samples with an even number of sites have shown that 
the ground state is a singlet and the existence of an exponentially low- 
lying triplet state (the so-called Kennedy triplet) in the Haldane gap. This 
leads to a fourfold ground-state degeneracy in the thermodynamic limit 
which can be simply understood by means of the VBS description of Fig. 
01 as described above. The existence of these two 5=1/2 edge states is 
a robust and striking property of the Haldane spin liquid phase and they 
are pres ent in the whole region < 1 of the bilinear-biquadratic spin- 
1 model The physical properties of these 5=1/2 chain-boundary 

excitatio ns iii the open spm-l chai n have been invest i gat ed in d etail nu- 
mericallyl ^^ l ^^ l '^^ l '^^ l '^^ l '^^ l '^^ l ^^ l ^ ^ and also analytically 1 ^ ^ | 82 |83] j^gniark- 
ably enough, this liberation of fractional spin-1/2 degrees of freedom in 
a spin-1 Heisenberg chain has been observed experimentally in the spin-1 
compound NENP cut by non-magnetic impurities (Zn^+ or Mg^+) l ^^ l g^v l 
and also doped with magnetic ions (Cu^"'") In particular, these spin- 
1/2 chain-boundary excitations reveal itself as satellite peaks in the NMR 
profile of the Mg-doped Y2BaNi05'^. A very recent inelastic neutron scat- 
tering experiments further probe the microscopic structure of these edge 
states through the wave-vector dependence of the Zeeman resonance. 



2.4. General spin-S case 

For higher spin-5 value {S > 1), the difference between integer-spin and 
half-integer AF Heisenberg chains has been explained within the semiclas- 
sical analysis pionne red by Haldane "2^*^ and also by means of the Abelian 
Eland non-Abelianllribosonization approaches. The Haldane's conjee- 
ture for S > 1 has also been investigated numerically and experimentally 
over the years. _______ 

The numerical studies I 91 l 92 | 9;j l 94 | ^^le S = 3/2 AF Heisenberg 
chain report that the model diplays a critical behavior that belongs to 
the same universality class as the S — 1/2 Heisenberg chain. Several 
quasi- I D AF compounds with spin 3/2 have been found like CsVBrs, 
CsVClal^^I^, and AgCrP2S6'^. In particular, the uniform susceptibility 
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of the spin-3/2 AgCrP2S6 material measured experimentally 1^ is consis- 
tent with the result for spin-1/2, i.e. confirming the universal behavior of 
half-integer spins. In the S = 2 case, the existence of a Halda ne g ap has 
been established in various qua ntum Monte-Carlo approaches and in 
DMRG calculations In the thermodynamic limit, a finite gap 

of A = 0.085(5) J and a fini te co rrelation length ^ = 49(1) lattice spacings 

have been found in DMRG EHI Several Haldane spin-2 compounds have 

1 9 

been discovered over the years (see Ref. for a review ) like CsCrCls and 
(2,2 -bipvridine)trichloromanganase(III) (MnCl3(bipy)) EH. Experiments 
| 1Uz |qj^ -j-j-^jg jg^g^ compound report the existence of a gap A/ J = 0.07 ±0.02 
which is in good agreement with the numerical results. 

Integer S > I spin chains share similar properties as the S = 1 case. 
In particular, there is still a hidden AF order in the integer S > 1 case. A 
non-local string order pararneter that measures this topological order has 
been proposed by Oshikawa^^ 



string 



I) = - lim (5f exp I itt/S \ Sf), (14) 




where a = x,y, z. This order parameter can be exactly computed in spin-^* 
gener alization of the VBS state I I I and numerically for the Heisen berg 
chain U ^ l ^^ l ^^^ l For instance, in the spin-2 case, DMRG calculations 
found O^^.j.jjjg (7r/2) = 0.726(2). However, the full characterization of the 
hidden AF order for integer spin S* > 1 is less clear than for S = \ where, 
as seen above, it is related to a broken Z2 x Z2 symmetry. In fact, the 
nature of the hidden symmetry breaking is still an open problem for S* > 1 . 
Insight on this issue might be gained by studying the structure of edge 
states of open integer spin-S chains. Higher-5' VBS states with open BC 
have (S -f- l)^-fold degenerate ground states with the presence of two free 
Sjl degrees of freedom at each end of the chain in close parallel to the S' = 1 
case (see Fig.OJ. A non- linear sigma model approach 'l^^ of the open spin-S* 
Heisenberg chain reveals the existence of spin-5'/2 edge states in the integer 
spin case. These results lead us to expect that Z2S x Z2S might be a good 
candidate to describe the hidden symmetry breaking scheme in the integer 
spin-S" case. In the 5 = 2 case, the approximate nine-fold degener acy of the 
ground state of finite open chains have been observed numericallyl^^l^^E^ 
by DMRG. The presence of these spin-1 chain-end excitations has also been 
confirmed experimentally in the ESR study of t he sp in-2 Haldane compound 
CsCrCls doped with non-magnetic Mg^"*" ions'''>i"J. 
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One- Dimensional Quantum Spin Liquids 15 

The liberation of these fractional spin-S'/2 degrees of freedom when non- 
magnetic impurities are introduced in integer spin-S* AF Heisenberg chain is 
a remarkable fact. It stems from the presence of a non-trivial topological or- 
dering in the ground state of the spin-S' chain. In this respect, the existence 
of these edge states opens the possibility to make a distinction between 
spin-S" Haldane liquid phases. These massive spin liquid phases share the 
same thermodynamics properties but have for instance different chain-end 
ESR response due to the special nature of their edge states. In fact, this 
topological distinction between gapped spin liquid phases can also be ap- 
plied to the massless case. All half-integer AF Heisenberg spin chains have 
the same bulk universal properties described by an IR fixed point which 
belongs to the su(2)i WZNW universality class. However, the ground state 
of half-integer spin chains still sustains a topologi cal o rdering. Indeed, by 
means of a non-linear sigma model approach, Ng CHSI has given evidence 
that edge states with fractionalized quantum number {S — l/2)/2 exist in 
open half-integer spin-S* AF Heisenberg chain. This interesting result has 
been confirmed numeri cally by examining open spin chains with s pins up 
to S = 9/2 by DMRG EJCnHi. Recently, a numerical investigation HHi of 
the non-local string order parameter 114(1 directly reveals the presence of 
a topological order in half-integer spin-S AF Heisenberg chains. The exis- 
tence of this hidden topological order (or the related chain-end excitations) 
leads to a topological distinction between gapless spin liquid phases with 
the same emerging quantum criticality. 



2.5. Two-leg spin ladder 

A further impetus for the study of low-dimensional spii i sys tems was given 
by the discovery of spin- ladder materials one decade ago^ '^ 1 These systems 
consist of a finite number (niog) of spin- 1/2 AF Heisenberg chains coupled 
by a transverse exchange interaction J^. The Hamiltonian of this model 
reads as follows 

ni(jg ni(jg — 1 

Ti-m^^ — J\\''^^''^^ Sa.i ■ Sa,i+1 + J±''^^ E/ ^a,i ■ Sa+IA, (15) 
i a—1 i a—1 

where Sa,i is a spin-1/2 operator at site i on the a-th chain and the in- 
chain exchange J\\ is antiferromagnetic (J|| > 0). One of the main interest 
of these spin ladders stems from the fact that they are intermediate objects 
between ID and 2D systems. Moreover, these quasi-lD systems share many 
properties with the cuprates and, because of their relative simplicity, may 
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provide some insights on mechanisms behind high-T^ superconductivity. 

Despite this original motivation to study spin ladders, these systems 
display striking properties which make them interesting by themselves in 
the field of ID quantum magnetism. In particular, in close parallel to the 
qualitative difference between integer and half-integer AF Heisenberg spin 
chains, the universal properties of spin ladders H15|l with open boundary 
conditions in the tr ansverse direction st rongly depend on the parity of the 
number of legs rtiog | lll | ll^ | llt^ | 114 | 1 15] Ladders with even number of legs 
are spin liquids with a finite gap in the excitation spectrum and expo- 
nentially decaying spin-spin correlations. The original massless spinous of 
the spin-1/2 chains get confined by the interchain coupling to form coher- 
ent optical spin-1 magnon excitations. Upon doping, these ladders display 
quasi-long-range s uper conducting pairing correlations with an approximate 
d-wave symmetry! " I In contrast, odd- legged ladders have a gapless spec- 
trum which is characterized by a central charge c = 1 corresponding to 
an effective S — 1/2 AF Heisenberg spin chain. Upon doping, odd-l egged 
ladders exhibit a metallic behavior typical of a Luttinger liquid This 
difference on the energy spectrum of the spin ladders depending on the 
parity of niog is strongly reminiscent of the Haldane conjecture for AF 
Heisenberg spin chains with S — niog/2. It can be simply understood from 
the fact that in the strong ferromagnetic rung coupling limit the spin ladder 
(|15|l is indeed equivalent to a single S — rticg/2 AF Heisenberg chain. 

Several spin ladders materials have been discovered over the years. The 
family of compounds Sr„_i Cun. j-i O^ tv. approximately realize ladders with 
"-leg = {n+\)/2, n = 3, 5, 7.. legs^ I The even-odd scenario in ladders was 
confirmed experimentally from susceptibility and muons r esonance mea- 
surements on this compound with 7i = 3 and n = 5 xhe existence 
of a spin gap in the two-leg case (riiog ~ 2) was found in experimen ts on 
the two-leg compo unds Ai 4 Cu^4 O4 1 (A14 = LafiC as, (Sr,La,Ca)i4) 
Cu2(C5Hi2N2)2Cl4^^^E2n (C5Hi2N)2CuBr4 and the organic com- 
pound BIP-BNO I '^ 1 However, it is worth noting that additional inter- 
actions are required to fully describe these materials such as for instance 
a non-negligeable ring exchange for the compound A14CU24O41 |1'^4 | 125 | 
or a small diagonal interchain exchange interaction for Cu2fC5Hi2N2)2Cl4 

The existence of a spin gap in the two-leg case can be anticipated by 
considering a strong-coupling analysis of the model H15|) with nicg = 2. 
In the limit of strong AF interchain coupling Jj^ 3> J\\ , the ground state 
consists of singlet bonds formed across the rungs of the ladder, with triplet 
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One- Dimensional Quantum Spin Liquids 17 

excitations separated by a large energy gap of order J±. Perturbations in 
small J|| will cause these singlet bonds to resonate and the triplet excitations 
form a band with bandwidth J|| but the spectral gap survives EEI in 
the strong ferromagnetic interchain —J± ^ J|| coupling case, local spins 
S = 1 associated with each rung of the ladder are formed leading thus to 
a spin-1 AF Heisenberg chain with a non-zero Haldane gap in the energy 
spectrum. The cross-over between strong and weak coupling li mits has been 
carefully analysed in numerical calculations | f ^7 | 128 | 129 | 130 | turns out 
that the ground state in the two strong-coupling limits | J±|/J|| » 1 evolves 
adiabatically with increasing J|| . The spin gap survives for arbitrarily large 
J||/|J^| and finally vanishes for J± = 0, i. e. when the two spin-1/2 AF 
Heisenberg chains are decoupled. The critical point J± — separates thus 
two strong-coupling massive phases: a rung-singlet phase for J± > and a 
phase with J± < which is smoothly connected to the Haldane phase of 
the spin-1 chain. 

The opening of the spin gap upon switching on the interchain coupling 
can be investigated by means of the bosonization approach -"^ In 
particular, the field theory that accounts for the massless spinon-optical 
magnon transmutation, when J± is small, corresponds to an SO (3) x Z2 
symmetric model of four massive Majorana (real) fermions i.e . four 

off-critical 2D Ising models (for a review see for instance the book^^. The 
resultin g low-en erKv field theory is described by the following Hamiltonian 
density 



- ^ {eBM% - eLd.a) - ^^se^fL + ^^arg, (16) 

where ^_r,l are a triplet of right and left-moving Majorana fermions that 
describe the 5=1 low-lying excitations of the two-leg spin ladder and the 
Majorana fermion accounts for the singlet excitation. The masses of 
these real fermions are given by in the weak coupling limit: mt — Jj_A^/27r 
and iris = — 3Jj_A^/27r, A being a non-universal constant. These triplet 
and singlet massive Majorana fermions are weakly coupled by a marginal 
perturbation associated to the last term of Eq. (|16(l : 

Hm^rs = ^ (^B. ■ ' + 92 6? ' C^Ci, (17) 



2 

with gi = —52 — 7raoJx/2, uq being the lattice spacing. The model Hlfi|l. 
derived in the weak coupling limit |J_l| <^ J||, is expected to have a larger 
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extent and to capture in fact the low-energy properties of the two-leg spin 
ladder for arbitrary Jx with a suitable redefinition of the masses mg^t, ve- 
locities Vs^ti and coupling constants gi 2- This stems from the continuity 
between the weak- and strong-coupling limits in the two-leg spin ladder 
observed numerically ^2^*12^*129 130 addition, in the strong ferromag- 
netic rung limit —J±^ J||, the singlet excitation described by the Majo- 
rana fermion ^ are frozen (|tos| — ^ oo) so that the low-energy properties 
of the model 1)16(1 are governed by the triplet magnetic excitations corre- 
sponding to the fermions ^_r,l- The resulting model coincides with the low- 
energy field theory of the spin-1 AF Heisenberg ch ain o btained by Tsvelik 
ESI by perturbing around the Babujian-Takhtajanll^ integrable point of 
the bilinear-biquadratic spin-1 chain (|5J) with /3 = — 1. The Hamiltonian 
(I16|l correctly thus captures the low-energy properties of the two-leg spin 
ladder in the limit —J± ^ J|| where it reduces to the spin-1 AF Heisenberg 
chain. 

Despite of its apparent simplicity, the model H16() is not an integrable 
field theory and takes the form a massive S0(3) x Z2 Gross-Neveu model 
However, the leading behavior of the physical quantities of the two-leg 
spin ladder can be determined b y tre ating the marginal contribution (|17() 
pcrturbatively. As shown in Ref. ' ^'^A ^ this term leads to a renormalization 
of the masses and velocities so that the low-energy description ((16|l simply 
reduces to four independent massive Majorana fermions, i.e. four decoupled 
off-critical 2D Ising models. This mapping onto off-critical 2D Ising models 
can then be exploited to derive the low-energy properties of the two-leg spin 
ladder such as, for instance, the spectrum of elementary excitations, low-T 
thermodynamics and the leading asymptotics of spin-spin correlation func- 
tions. To this end, one needs to express the spin operators Sa,i in terms of 
the Ising fields. In the continuum limit, the spin densities separate into the 
smooth and staggered parts'^ Sa{x) — JaL{x) + JaR{x) + {—lY^°'°na{x). 
The chiral su(2)i currents, uniform parts of th e spin de nsities, can be writ- 
ten locally in terms of the Majorana fermions I^^ESl 

JaR,L = - ^ A + «y ^R,lS,R,L: (18) 

with Ti = +1 and T2 — —1. In contrast, the staggered magnetizations fia 
are non-local in terms of the underlying fermions a nd express i n terms of 
the order (ctq) and disorder (^a) of the Ising models I^SEIESI 

n+ ^ i^J.a^ila■2cr3, ^J,oO■l^i20■3, ^J■acrlCr2^J■3) 

ft- ~ CT0/XlCr2/X3, Cro/Xl/Z2Cr3) , (19) 
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One- Dimensional Quantum Spin Liquids 19 

where n± = ni ± n2. Due to the existence of non-zero masses rris^t for the 
fermions in Eq. (|16(l . all related Ising models are non-critical. A spectral 
gap is thus present in the two-leg spin ladder for all signs of the interchain 
interaction Jx and it opens linearly with J± in the weak coupling limit. 
In addition, the Ising description H18I19() of the spin densities allows the 
calculation of the leading asymptotics of spin-spin correlations using exact 
results of the two-point function of a non-critical Ising model. Since the signs 
of the tripl et and singlet masses are always opposite in Eq. IjlGII . it can be 
shown I I that the dynamical spin susceptibility displays a sharp single- 
magnon peak near q — 7r/ao and lo — |mf| which reflects the elementary 
nature of the triplet excitations for all signs of J±. 

In this respect, the spin liquid phase of the two-leg spin ladder is very 
similar to the Haldane phase of the spin-1 chain. However, it has been 
stressed recently that the phase with J± < 0, smoothly connected to the 
Haldane phase of the spin-1 chain, and the rung-singlet phase for Jx > 
are, in fact, topologically distinct gapped phases -121.. The distinction is in- 
timately related to the short-range valence bond structure of their ground 
states. In close parallel to the classification of short -range valence bond con- 
figurations on a two-dimensional square lattice^^ two different topological 
classes can be defined in the one-dimensional case by c qunti ng the number 
Qy of valence bonds crossing an arbitrary vertical line l ^^^ l In the case of 
the rung-singlet phase, Qy is always even while it is odd for the Haldane 
phase with J± < 0. Two different non-local string order p arameters ca n 
then be defined in connection to this topological distinction EHEO^SZl 

= - hm ( {Sl,+, + S^,,) S r4+i (si.+.+si, ) ^ga ^ ^ 

\t~j\—>--\-oo 

a"dd = -, lini ((5M + 52".0e*"^™(^°'=+^^°^)(5r.,+52",.^ (20) 

\t—j\—>--]-oc 

with a ^ x,y^ z. The Haldane and rung-singlet phases are then character- 
ized by OZ^^ + 0, = 0, and by ^ 0, O^^d = respectively. The 

order parameters 12()|l reveal also the difference nature of the hidden AF 
Neel order in the Haldane phase (non-zero triplet states along the rung) 
and in the rung-singlet phase (non-zero triplet states along the diagonal) 
[l29 130 137* This topological difference can be discussed in light of the Ising 
description (|19|l of the two-leg spin ladder. The string order parameters 
can be expressed in terms of the order and disorder Ising operators ' 
©odd (o'i)^(o'2)^ and Ocvon ~ (/^i)^(M2)^- For a ferromagnetic (respec- 
tively antiferromagnetic) interchain coupling, the Ising models in the triplet 
sector are in their ordered (respectively disordered) phases so that Oo^d 7^ 
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and Oevon = (respectively O^^^ = and C|ven 7^ 0). In this Ising descrip- 
tion, the phases with J± > and J± < are thus simply related by a 
Kramers- Wannier duality transformation on the underlying Ising models. 

The topological distinction between these two phases becomes manifest 
when analysing the ground-state degeneracy depending on t he n ature of 
BC used. In the open BC noted by the authors of Ref . ground 

states of gapped spin liquid states characterized by an odd value of Qy have 
spin-1/2 edge states, while these end states disappear when Qy is even. The 
existence of these S = 1/2 chain-end degrees of freedom leads to a ground- 
state degeneracy in a two-leg spin ladder with open BC. In the Haldane 
phase of the ladder with J± < 0, finite open chains with an even number 
of sites have a singlet ground-state with an exponentially low-lying triplet 
in the spin gap resulting on a fourfold ground-state degeneracy in the ther- 
modynamic limit. In contrast, the ground state in the rung-singlet phase is 
always unique whether open or periodic BC are used: no low-lying triplet 
states are found inside the spin gap in open ladder with J± > and an 
even number of sites. In this respect, the rung-singlet phase with J± > 0, 
in contrast to the J± < case, is not equivalent to the Haldane phase char- 
acterized by 5 = 1/2 chain-end excitations even though they share similar 
properties such as the presence of a spin gap, and a non-zero string order 
parameter. The existence or absence of spin-1/2 edge states in the open two- 
leg spin ladder can also be discussed within the Ising model description (fT?)|l 

In particular, the fourfold ground-state degeneracy of the open two-leg 
ladder with a ferromagnetic interchain coupling can be obtained within this 
approach. Indeed, for J± < 0, the three Ising models for the triplet sector 
are all in their ordered phases while the Ising model for the singlet degrees 
of freedom belongs to its disorder phase so that (ct^) ^ (« = 1,2,3) and 
((To) = 0. In that case, each Ising model in the triplet sector has a doubly 
degenerate ground state which gives thus an eightfold degeneracy. However, 
there is a redundancy in the Ising description since the triplet Hamiltonian 
in Eq. H16|l . the total uniform and staggered magnetizations in Eqs. (|18I19(I 
are all invariant under the transformation: Cr l ^ "^RLTt^i ~* M'j ^^'^ 
(Tj — CTi, i = 1,2, 3. This leads to a physical fourfold ground-state degen- 
eracy as it should be in the Haldane phase. Another advantage of this Ising 
model description is to make explicit the spontaneous breaking of a hidden 
Z2 X Z2 symmetry associated to the ground-state degeneracy for J± < 0. 
The existence of this hidden symmetry can also be revealed with help of a 
non- local unitary transformation on the l atti ce spins of the two-leg ladder 
EHIas in the spin-1 AF Heisenberg chain 
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Finally, the Haldane and rung-singlet phases of the two-leg spin ladder 
can also be distinguished by their optical properties. In the rung-singlet 
phase, the existence of singlet or triplet two-magnon bound states below 
the two-magnon continuum has been predicted theoretically by a number 
of groups 'i^^ 14 2 143, lH] ^ Recently, the singlet two-magnon bound state 
has been observed in t he optical conductivity spectrum of the compound 
(La,Ca)i4Cu2404i li^SI which contains layers with CU2O3 two-leg spin lad- 
der with Jj_ > 0. On the contrary, it is expected that for a ferromagnetic 
interchain coupling no bound states are presen t be low the two-particle con- 
tinuum as in the spin-1 AF Heisenberg chain ESI The existence of bound 
states in the two-leg spin ladder can be discussed in the context of the field 
theory approach 1)16(1 . Indeed, the marginal term (|17|l plays its trick by giv- 
ing rise to an effective interaction between the magnon excitations. In the 
AF interchain case, the contribution (|17|l is in fact marginal relevant and the 
interaction between magnons are attractive. The presence of bound states 
in this case can be argued qualitatively as a result of this attractive interac- 
tion. The Hamiltonian H17|l takes the form of a SO (4) Gross-Neveu model 
up to irrelevant contributions and a duality tran sfor mation —^ — in 
the singlet sector. The latter model is integrablel^ and, for J± > 0, the 
low-lying excitations are massive kinks and anti-kinks interpolating between 
the two degenerate ground states resulting from the spontaneous symmetry 
breaking of the Z2 symmetry: ■f^^f, "Cfl/^Li = 0, .., 3. The mass terms 
in the full Hamiltonian 1(16(1 break explicitely this Z2 symmetry and the 
ground-state degeneracy is lifted so that kink configurations are no longer 
asymptotic states of the field theory. The situation is in close parallel to the 
two-dimensional Ising model in its low-temperature phase upon switching 
on a magnetic field. The mass terms are then expected to induce a linear 
confining potential between the kinks giving rise to a sequence of bound 
states. In contrast, when J± < 0, the perturbation 1(17(1 is a marginal irrel- 
evant contribution and the spectrum of the model ((16() consists of massive 
fermions and their multiparticle excitations. 



2.6. Non-Haldane spin liquid 

The Haldane and rung-singlet phases of the two-leg spin ladder are not the 
only possible gapped spin-liquid states available in ID u nfrustrated quan- 
tum magnets. In this respect, Nersesyan and Tsvelik^^ have discussed the 
example of a gapped spin liquid phase without any coherent magnon exci- 
tations. The spectral function of this state displays a broad feature rather 
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than a single sharp magnon peak. Such a spin liquid can be stabilized by the 
introduction of a four-spin interchain interaction that couples two spin-1/2 
AF Heisenberg chains: 

2 

Hso = J ^ ^ Sa,i ■ SaA+1 + K ^ (^SlJ ' Sl^i+lj ^5*2,; • 5*2, . (21) 
i a—1 i 

The biquadratic interaction represents an interchain coupling for the spin- 
dimerization operators {ea,i ^ (— l)*'S'a,i+i ■ Sa^i) of each chain which can be 
effectively generated by spin-phonon interaction. A second motivation to 
investigate the effect of this four-spin interaction stems from orbital degen- 
eracy. In most of transition metal compounds, in addition to the usual spin 
degeneracy, the low-lying electron states are also characterized by orbital 
degeneracy A starting point to study magnetic properties of magnetic 
insulators with Jahn- Teller ions is the two-band Hubbard-like models. At 
quater- filling (one electron per atom), this system is a Mott insulator in the 
limit of strong Coulomb interaction and the state of each ion can be char- 
acterized by a spin degrees of freedom Si^i and an orbital state described 
by a pseudo-spin- 1/2 52, i. In the large Coulomb repulsion limit, the sim- 
plest Hamiltonian that describes the competition between spin and orbital 
degrees of freedom in one dimension reduces to the spin-orbi tal m odel H21|l . 
In particular, this model has been introduced by Pati et al. ^^to explain 
the unsual magnetic prop eties of the quasi-one-dimensional spin gapped 
material Na2Ti2Sb20'2Hl. 

The spin-orbital model H21|l is unfrustrated for i^T < whereas the 
frustration manifests itself in the antiferromagnetic case K > only in the 
intermediate regime iiT ~ J as we shall see later. The Hamiltonian (|21|) is 
invariant under independent SU(2) rotations in the spin (Si) and orbital 
(^2) spaces. For generic couplings, the model (|^ is thus SU(2) x SU(2) x 
Z2 symmetric, the additional Ising symmetry being the exchange between 
the spins 5*1 and ^2. In the weak coupling limit \K\ <C J, this underlying 
SO (4) symmetry of the spin-orbital model (|21|l is reflected in the form of its 
low-energy Hamilt onia.n den sity which is described in terms of four massive 
Majorana fermions | 1^^ | 149 | 

Hso =^ ~y E (^R^-^R - -^^T. ^R^i' (22) 

a=0 a=0 

with m = aK, a being a positive non- universal constant. In Eq. H22|l . a 
marginal contribution, similar to Eq. (|17|l . has been neglected. In contrast 
to the low-energy description l(T()| of the standard two-leg spin ladder, the 
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triplet and singlet Majorana modes become equally important here as a 
consequence of the SO (4) symmetry of the spin-orbital model (|21|l . The re- 
sulting low-energy properties of the model can then be determined by this 
mapping onto four off-critical 2D Ising models. As the staggered magneti- 
zations (|19ll . the spin dimerization fields ei.2 can be expressed in terms of 
the order and disorder operators of the underlying Ising models: 

Ca ~ ^1^2/^3/^4 ± Ta<7ia2(J3a4- (23) 

For an AF biquadratic coupling {K > 0), the four Ising models belong to 
their ordered phases so that the model H22|l enters a spontaneously dimer- 
ized phase with a finite gap and (ei) = —{^2) = i|eo|- The ground state 
is thus two-fold degenerate and dimerizes with an alternating pattern as 
shown in Fig. 01 (a). In contrast, in the ferromagnetic K < case, the 

(a) 

» ^ » * » # 



3 » i> » * £ 



(b) 

(» m m B m i> 



(i m m "» m T> 

Fig. 4. Ground states of the spin-orbital model in the weak coupling limit \K\ <C J- 
The bonds indicate a singlet pairing between the spins, (a) Alternating dimerization 
for K > 0; (b) in-phase dimerization for ii" < 0. In each second ground state 

is obtained by applying the one-step translation symmetry on each chain to the states 
depicted here. 

dimerization is now in-phase between the two chains: (ei) — (£2) — ±|eo| 
as depicted by Fig. 01 (b). A first distinction between this gapped spin 
liquid state and the Haldane and rung-singlet phases of the two-leg spin 
ladder stems from the ground-state degeneracy and the fact that the lat- 
tice translation symmetry is spontaneously broken in the ground states 
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of Fig. ^ A more drastic difference appears at the level of the low-lying 
excitations. Indeed, the elementary excitations of the spin-orbital model 




Fig. 5. Triplet excitation of the alternating dimcrization described in terms of a pair of 
dimorization kinks of each chain. 

are neither optical magnons nor massive spinons but a pair of propagat- 
ing massive tr iplet or singlet kinks connecting two spontaneously dimerized 
ground states 'J'^. An example of such a triplet excitation is described in 
Fig. O for the staggered dimcrization phase with K > 0. The composite 
nature of the S = I excitation reveals itself in the dynamical structure 
factor of the model which can be determined using the Ising description 
(|19|l of the total and relative staggered magnetizations-!-— . The dynamical 
magnetic susceptibility displays a two-particle threshold instead of a sharp 
magnon peak near q — tt/oo as in the two-leg spin ladder or the spin-1 
chain. This incoherent background in the dynamical structure factor leads 
to a new gapped spin liquid phase in unfrustrated quantum magnetism 
without any coherent magnon excitations. In this respect, this stat e has 
been called non-Haldane spin liquid by Nersesyan and Tsvelik ESI xhe 
existence of this spin liquid phase has been confirmed non-perturbatively 
at two special points in the phase diagram of the spin -orbital model H21|l . 
At J = 3K/4:, the ground state is exactly known '■'■^ and is a product of 
checkerboard-ordered spin and orbital singlets as in Fig. 01 (a) obtained in 
the weak coupling limit. At K = — 4J, the model is exactly solvable and 
the g roun d-state energy and triplet energy gap have been determined ex- 
actly 11^. The system has two spontaneously dimerized ground states and 
belongs to the class of non-Haldane spin liquid. 

Finally, it is worth noting that the alternating and in-phase dimcrization 
phases can be distinguished in close parallel to the topological distinction 
between the two gapped phases of the two-leg spin ladder. Using the topo- 
logical criterion of Kim et al. -^"^^ one observes from Fig.0]that the number 
Qy of valence bonds crossing an vertical line is odd (respectively even) for 
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the staggered (respectively in-phase) dimerization. The two non-Haldane 
spin hquid phases with K > and K < belong thus to two different 
topological classes. This difference becomes manifest when investigating 
the structure of the edge states of the two phases with open BC. The phys- 
ical properties of these boundary excitations of the semi- infinite spin-orbital 
model (|21|1 can be determined by means of the Ising mapping H22|l of the 
weak coupling limit similarily to the cut two-leg ladderl^. For K < 0, such 
edge states are absent w here as, in the K > case, two spin-1/2 chain-end 
excitations are expected H^^. 



3. Frustration effects 

In this section, we shall review some of the main aspects of the interplay 
between frustration and quantum fluctuations in AF spin chains and spin 
ladders. The natural question is whether frustration can stabilize new types 
of spin liquid phases with exotic spin excitations not encountered in Section 
2. The paradigmatic model to analyse the effect of frustration in spin chains 
is the Jl — J2 spin-1/2 Heisenberg chain with Hamiltonian: 

H^JlJ2s^■ S,+i + J2 5] 5, • S,+2, (24) 

i i 

where the next-nearest neighbor coupling J2 > is a competing AF in- 
teraction which introduces frustration. This model can also be viewed as a 
frustrated two-leg spin ladder where the spin chains are coupled in a zigzag 
way as shown in Fig.|Bl Let us first discuss some of the main characteristics 
of frustration in ID spin systems by means of a semiclassical approach. 




Fig. 6. Two-leg zigzag ladder. 



3.1. Semiclassical analysis 

Classically, one of the main effect of frustration is to favor non-coUinear 
magnetic ordering where the spins lie in a plane rather than along a single 



February 2, 2008 2:33 



WSPC/Trim Size: 9in x 6in for Review Volume 



revspinliquid 



26 

direction as in unfrustrated spin systems: 

Si ^ S cos (q ■ ni + S sin {q ■ n2, (25) 

where r; is the spatial location of the site i and ni^2 are two mutually orthog- 
onal fixed vectors in spin space with unit length. The magnetic structure 
corresponding to the ordering H25|l is a circular spiral with a pitch angle 
related to the wave-vector Q which is, in general, incommensurate. In the 
case of the Ji — J2 chain 1)24(1 . the spiral ground state is stabilized when 
J2/J1 > 1/4. The spins are arranged in a canted configuration in which 
each spin makes an angle a with its predecessor such that cos a = — Ji/4 J2. 
The classical ground state of the model is doubly degenerate since the spin 
configurations (|25|l can turn clockwise and counterclockwise with the same 
energy along the spiral axis 713 = ni A r?2 . A Z2 discrete symmetry charac- 
terized by this right- and left-handed chirality is thus sp ontan eously broken 
in this helical structure. A corresponding chiral ordering^SI can be defined 
and detected by the chiral order parameter: {{Si A S'i+i) • ^13) =/= 0. In con- 
trast to the Neel state of unfrustrated magnets, the SU(2) spin rotation of 
the Heisenberg Hamiltonian is completely broken in the non-collinear state 
()25|l . Three Goldstone or spin- wave modes are thus expected here from this 
spontaneous symmetry breaking scheme. It leads to a CFT with central 
charge c — S with extended criticality in comparison to the unfrustrated 
case (c = 2jSee Section 2). A second distinction stems from the spin- wave 
calculation'-"^ for non-collinear ordered states which reports the presence 
of two different spin-wave velocities. Therefore, the low-energy field theory 
describing canted magnets is expected to be non-Lorentz invariant. To iden- 
tify this effective field theory, one needs to define a suitable order parameter 
for canted magnets. The long-wavelength fluctuations of the spin configu- 
rations ((25() can be captured by introducing a two-component complex field 
Za,a = 1,2 of unit modulus (l^ip -I- |z2p = 1) such that^l 

ni + OT2 = eacZc^abZb, (26) 

where Sab is the standard antisymmetric tensor and is a vector formed by 
the Pauli matrices. It is indeed straightforward to see that this identification 
((26|l correctly reproduces the constraints: n\ = n^ = l and fii-n2 = if the 
two-component complex field Za belongs to the three-dimensional surface 
of the unit sphere in four dimensions denoted by S'^. Moreover, it is worth 
noting that the field z = (21,2:2) transforms like a S = 1/2 spinor under 
spin rotations. The representation is double-valued since z and —z 
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describe the same non-collinearly ordered state. The order parameter space 
corresponding to canted magnets H25(l is thus 

A crucial consequence of this identification is that the order param- 
eter allows topologically nontrivial vortices having a 2it circulation since 
TTi (S ^ /Zy) = ^2'!^. Upon encircling such a Z2 vortex, called now a vison 
by a closed loop, the value of the spinor complex field z = (21,2:2) 
changes smoothly from z to — z, i.e. two points diametrally opposed on the 
surface of a sphere S'^ . The existence of these stable vison defects reveals it- 
self in the semiclassical description of frustrated spin chains like the Ji — J2 
Heisenberg chain H24|l . In particular, in sharp contrast to unfrustrated spin 
chains, no Pontryagin topological term has been found in the effective ac- 
tion that desc ribes the low-energy properties of the Ji — J2 Heisenberg 
ehainHSmSHl 

and two-dimensional frustrated Heisenberg antiferromagnets 
| 159 |^j^ ^Y^Q large spin limit. This stems from the very special topological na- 
ture of the order parameter for canted magnets which has a trivial second 
homotopy group: 'K2{S^ /Z2) — ^"^^ . However, there is still a subtle distinc- 
tion bet ween half-integer and integer spins in the Ji — J2 Heisenberg chain 
[160 157 r ela-ted to the existence of stable visons. A Berry phase calculation 
[I 57 100 l(jl | l(?2 | gijows that tunneling between sectors with different Z2 to- 
pogical number is possible for integer spins but not for half-integer spins 
due to cancellation between pairs of paths. The ground state of the Ji — J2 
model is thus non-degenerate for integer spins and two-fold degenerate (two 
different sectors non-coupled by tunneling effect) for half-integer spins. The 
quantum number carried by el ementary excitations depends also strongly 
on the nature of the spins'""^. In the half-integer case, there is an energy 
gap towards the creation of Z2 visons so that these topological defects are 
strongly suppressed. The low-energy excitations are then described by the 
spinor field z — (zi, 22) of Eq. (|26(l which is free to propagate in absence 
of visons and carries a 5* = 1/2 quantum number under spin rotations. Ne- 
glecting anisotropy, the Euclidean effective action that governs the leading 
low-energy properties of the Ji — J2 chain for half-integer spins is the SU(2) 
X SU(2) principal chiral model defined by 

5g = J y dxdr Tr {d^G^d^G) , (27) 

where G is a SIJ(2) matrix formed by the two complex fields 21.2. This 
field theory l l^^^l^^^ l integrable and its low-energy spectrum consists of 
massive excitation with 5* = 1/2 quantum number i.e. massive deconfined 
spinous. In contrast, in the integer spin case, the visons now proliferate in 
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the ground state so that the spinor field z = {zi, Z2) cannot any longer be 
defined as single- valued configurations. The suitable fields that describe the 
low-energy excitations of the Ji — J2 chain with integer spins are bilinear 
of z and correspond to the ni_2 fields of Eq. if^ . In that case, neglecting 
anisotropy, the effective field theory that captures the leading low-energy 
properties of the frustrated spin chain with integer spins is the non-linear 
sigma model SO (3) x SO (3) with action: 



where i? is a rotation matrix made of the triplet vectors ni^2 and 713 = 
ni A ^2 • The low- l ying spec trum of this field theory can be determin ed by 



it consists of a triplet of = 1 massive magnons. In the integer spin case, 
the spinous are thus expected to be confined into optical magnons. 

An important prediction, obtained in this semiclassical description, is 
thus the existence of a gapped spin liquid phase for half-integer spins with 
a two-fold degenerate ground state and massive deconfined spinons. Such a 
state is stabilized by frustration and represents a spin liquid phase not en- 
countered in Section 2. We shall now consider the ultra-quantum case with 
S" = 1/2 where a combination of numerical and field theoretical techniques 
can be used to fully determine the main characteristics of this spin liquid 
phase. 

3.2. Spin liquid phase with massive deconfined spinons 

The phase diagram of the spin-1/2 Ji — J2 Heisenberg chain H24|l has been 
studied extensively over the years after the bosonization analysis of Haldane 
. This problem is not als o a pure ly academi c que stion since inorganic 
compounds such as CuGeOa l ^vP l lPi l LiV205 l -^vffl can be considered as 
prototypes of the spin-1/2 Ji — J2 chain H24|l . In particular, values such 
as Ji ~ 160 K and J2/J1 — 0.36 have been proposed for CuGeOa 
In addition, the quasi-lD compound SrCu02 contains a collection of spin- 
1/2 Heisenber g chains assembled pairwise in an array of weakly coupled 
zigzag 

ladders These zigzag ladders are built from corner-sharing 

Cu-0 chains with an exchange J2 staked pairwise in edge-sharing geome- 
try. The frustrating interaction Ji between the chains (see Fig. stems 
from the nearly 90° Cu-O-Cu bonds and is expected to be weak 11^. Fi- 
nally, two other possible realizations of the spin-1/2 Ji — J2 chain (|24|l 
have been recently proposed: the compound (N2H5)CuCl3 which can be 




(28) 



a large N approach 




or by means of a strong-coupling analysis 
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described as a two-leg zigzag spin ladder with J1/J2 ~ 0.25^21 and Cu[2- 
(2-aniinomethyl)pyridine]Br2 with a ratio J2/J1 — 0.2^^. 

A starting point for investigating the phase diagram of the spin- 1/2 
Ji — J2 chain H24|) is to consider the weak coupling limit when J2 ^ Ji. 
This enables us to study the stability of the massless spinous of the spin- 
1/2 AF Heisenberg chain upon switching on a small next-nearest-neighbor 
frustrating interaction J2. In this weak-coupling regime, the low-energy 
effective Hamiltonian density of the model reads as follows ^3 

Hcff = ^ ( + 4) + jJl ■ Jr, (29) 

where Jl,r are the left and right su(2)i currents that generate the su(2)i 
quantum criticality of the spin-1/2 AF Heisenberg chain. The low-energy 
physics of the model are thus mainly determined by the marginal current- 
current interaction of Eq. (|29|) with coupling constant 7 ~ J2 — J2C, J2c 
being a non- universal positive constant. For a small value of J2, one has 
7 < so that the interaction in Eq. 129|l is a marginal irrelevant con- 
tribution. The low-energy physics of the Ji — J2 chain is thus identical 
to that of the spin-1/2 AF Heisenberg chain which is governed by the 
su(2)i fixed point with additional logarithmic corrections introduced by 
the current-current interaction. In this respect, frustration plays no impor- 
tant role in the regime J2 < J2C and the spin-1/2 Heisenberg phase with 
massless spinons is stable upon switching on a small value of J2. However, 
for J2 > J2C (7 > 0), the current-current interaction becomes marginal rel- 
evant and a strong coupling regime develops with a dynamically generated 
spectral gap. A phase transition of Berezinskii-Kosterlitz-Thouless (BKT) 
type Cil occurs at J2 = J2C which separates the gapless spin-1/2 Heisen- 
berg phase from a fully massive region. The actual value of the transition 
--P.^I^ ^^^'^ numerically by different groups to be: — 0.2411Ji 
| 175 | 17(j | 177 | 17§l At this point, it is worth noting that the quasi-lD material 
Cu[2-(2-aminomethyl)pyridine]Br2 can be described by a spin-1/2 Ji — J2 
Heisenberg chain with J2/J1 — 0.2 and should thus belong to the crit- 
ical Heisenberg phase with J2 < J2c- The absence of a spin gap for this 
compound has been reported experimentally . 

The main characteristics of the strong coupling massive phase with 
J2 > J2C can be determined from Eq. (|29(l . The field theory H29|l is indeed 
integrable and corresponds to a chiral Gross-Neveu model or a non-Abelian 
version of the Thirring model with a SU(2) symmetry. The low-energy ex- 
citations of the model (^i^ are a massive doublet (massive spinons) with 
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mass m 



a. 



'0 



exp 

(_27r/7)iniiHni. 

A striking effect of frustration is thus 



the formation of a mass for the spinous of the spin-1/2 AF Heisenberg chain 
without confining them into S* = 1 excitations as in the two-leg spin lad- 
der. A second consequence of frustration is the presence of spontaneously 
dimerization in the model. A simple way to exhibit this dimerization is to 
use the bosonization approac h to e xpress the Hamiltonian (|29|l in terms of 



where $_r,l are the chiral components of the bosonic field $ ($ = ^b + ^l) 
and Q is its dual field {Q = — ^r)- The SU(2) invariance of the model 
()30|l is hidden in the structure of the interaction with a single coupling 
constant and the fact that the bosonic field $ is compactified on a circle 
with a special radius R — l/\/27r consistent with the SU(2) symmetry. This 
compactification leads to the following identification: 



n being integer. In the phase with J2 > J2C, one has 7 > so that the 
bosonic field is pinned at one of its minima: ($) = py^ n/2, p being integer. 
However, by taking into account the identification l|31() . the = Stt sine- 
Gordon model (|30|l has only two inequivalent ground states with ($) ~ 
and ($) = \/tt/2. This two-fold degeneracy can be interpreted as resulting 
from the spontaneous breaking of a discrete Z2 symmetry. This symmetry 
identifies with the one-step translation s ymm etry (Tag) which is described 
by the following shift on the bosonic field + \/tt/2. This symme- 
try is spontaneously broken in the phase J2 > J2C and the two ground-state 
field configurations are connected by this translation symmetry. An order 
parameter designed to characterize this phase is the spin dimerization oper- 
ator e„ = (— • ^n+i which admits the following bosonic representation 
in the continuum limit e ~ cos-\/27r$ ■^^ This operator changes sign un- 
der the lattice translation symmetry and has a non-zero expectation value 
(e) 7^ in the two ground states of the = Stt sine-Gordon model Ip^Hjl . 

For J2 > J2C, frustration stabilizes thus a gapful spontaneously dimer- 
ized phase which is characterized by a two-fold degenerate ground state 
and a spontaneous breaking of the lattice translation symmetry. The ele- 
mentary excitations of this phase are massive spinous that carry S* = 1/2 
quantum number and identify with the kinks of the underlying dimeriza- 
tion. A simple way to understand the emergence of these deconfined spinous 



a = Stt sine-Gordon model ^5 



cos V87r$, (30) 




(31) 
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is to consider the Majumdar-Ghosh (MG) point at J2 = Ji/2 where 
the ground state of the lattice model (|24|l is exactly known. This MG point 
for the spontaneously dimerized phase plays a similar role than the AKLT 
point (|10|) for describing the main properties of the Haldane phase in the 
phase diagram of the spin-1 bilinear-biquadratic chain Q. For J2 — Ji/2, 
the Hamiltonian (|24|) takes the following form up to a constant: 

O T 

^MG = ^^^'3/2(i-l,*,i + l), (32) 

i 

where P3/2 (i — + 1) projects the total spin of the three spins located 
at sites z — + 1 onto the S — 3/2 subspace. For an even number of 
sites N and periodic BC, the ground state of the MG model H32|) is two-fold 
degenerate and corresponds to the two singlet states: 

|a>i) = [1,2] [3, 4] ...[TV -1,7V] 

|a>2) = [2,3][4,5]...[7V,1], (33) 

where [i, z -I- 1] denotes a singlet bond between the spins at the i and i + 1 
sites. The two ground states are represented in Fig. [7| The lattice 



(a) 
(b) 



Fig. 7. The two dimerized ground states at the Majumdar-Ghosh point; (a) (b) 
l*2>. 



translation symmetry 7^^ is broken in the two ground states (|33|l and it 
exchanges them: 7^o|$i,2) = |'I'2,i)- The existence of gap in the MG model 
(|32|l has been shown rigorously by Affleck et al. and the spin- spin cor- 
relation function can be determined exactly at the MG point I t 

^ ^ 3 3 
{Si ■ Sj) = -Si^j - gi^ii-ji,!. (34) 
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This correlation function is thus zero for distance larger than one lattice 
spacing. The ground-state degeneracy and the spontaneous breaking of a 
discrete Z2 symmetry suggest the existence of topological excitations which 
interpolate between the two ground states These kink excitations can 
be viewed as the insertion of a spin- 1/2 in a sea of singlet valence bond 
states as shown in Fig.|Sl This spin-1/2 excitation (spinon) is nothing but 
a domain wall between the two states For periodic BC, these kinks 



, 1 



> » * € 



Fig. 8. Dimerization kink (spinon) with S = 1/2 quantum number. 

states appear always in pairs (see Fig. O and the low-l ying excitations of 
the MG model can be built starting from the state ^2 

|p,m) = [l,2]..[2p-3,2p-2] a2p-i[2p,2p + 1] ... 

[2m -2,2m- 1] aa-m [2m + 1, 2m -|- 2] ... [N - 1, iV] , (35) 

where a2p-i and a2m denote spin-1/2 states located at the sites 2p— land 
2m respectively. A variational approach of the low-lying excitations^^^SI 
can then be done by considering a linear combination of states (I35II and it 
leads to a triplet and singlet continuum spin excitations. In particular, 



» » 1 » < 1 » < 



Fig. 9. Spin-1 excitation built from two spinous. 

the dispersion relat ion of a massive spinon, obtained within this approach, 
takes the formCSl 

+ 4 cos(2fcao). (36) 

o Z 
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This dispersion relation has been verified numericaUy by exact diagonahza- 
tions Finally, as shown by Gaspers and Magnus ^^•-'^ there are addi- 
tional exact singlet and triplet bound-states at momentum k = 7r/2ao which 
are degenerate with energy E = Ji. Exact diagonalization and DMRG cal- 
culations ='-°'*'have confirmed the existence of these bound-states for a small 
range of momenta close to k = 7r/2ao in the region J2 > ^i/2. 

In summary, the spontaneously dimerized phase for J2 > J2C repre- 
sents a distinct spin liquid phase stablized by frustration. The main distinc- 
tion, from the spin liquid phases of the two-leg spin ladder and spin-orbital 
model, originates in the fractionalized nature of the quantum number car- 
ried by elementary excitations. In particular, instead of a sharp S = 1 
magnon peak as in the two-leg spin ladder, the dynamical structure factor 
of the model in this dimerized phase displays an incoherent background with 
additional bound-states features. The main difference between the dimer- 
ized phase of Fig.[7|and the staggered dimerization phase of the spin-orbital 
model (I21|l stems from the nature of the low-lying excitations. In the lat- 
ter phase, the excitations are a pair of propagating massive triplet kinks, 
as described above, whereas here the elementary excitations are massive 
spinous. 

This dimerized phase, with deconfined massive spinous excitations, ex- 
tends in the entire region with J2 > J2C as it has been shown numerically. 
Fig. El represents the evolution of the spin gap A, computed by DMRG 
as function of the next-nearest neighbor interaction J2. As depicted 




Fig. 10. Evolution of the spin gap, computed by DMRG, as a function of J2 for J2 > J2c', 
taken from Ref. 178. 
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by Fig. 1101 the spin gap is maximum for Jij J\ ~ 0.6 and decreases to zero 
in the large J 2 Umit. The absence of a spin gap in this regime can be easily 
understood since for J2 S> Ji , the J\ — J2 spin chain can be better viewed 
as a two-leg zigzag ladder (see Fig.EJ. In the limit J2 — > +00, the two spin- 
1/2 Heisenberg chains are decoupled so that the model is critical: A = 0. 
From the results of Fig. El it is very tempting to conclude that the whole 
region with J2 > J2C describes a single phase. However, a subtle qualitative 
change occurs in the model at the level of spin correlation functions after 
the MG point in close parallel to the Haldane phase ofthespin-l bilinear- 
biquadratic chain ® just after the AKLT point (0 '^31201113. Indeed, an 
incommensurate behavi or d evelops in the r eal-space spin-spin correlation 
function for J2 > Ji/2 [TT6,188,178|189|190 | j,^^^ incommensurability can 
be interpreted as a quantum signature of the spiral structure of the classical 
ground state of the Ji — J2 spin chain for J2 > Ji/4. After the MG point, 
the leading asymptotic of the spin-spin correlation function behaves as 

{Si ■ Si+r) ~ -f/j cos {q rao) exp (-rao/C) , (37) 

where the oscillation factor depends on J2/J1. In parallel to the classifi- 
cation of commensurate-incommensurate transitions in Ising spin systems 
the MG point, as the AKLT point is a disorder point of the 

first kind The momentum ^max which maximizes the static structure 
factor (Fourier transform of the spin-spin correlation function) become s in- 
commensurate not at the MG point but slightly after J2L — 0.5206Ji^SI 
at a Lifshitz point ' . For J2 > J2L , the static structure factor displays 
a double peak structure rather than a single peak. The low-lying excita- 
tions are then characterized by an incommensurate momenta gmax 7^ t^/o-o- 
In contrast to the classical case, it is worth noting that the BKT phase 
transition between the critical and dimerized phases and the onset of in- 
commensurability at the Lifshitz point occur at different points of the phase 
diagram, J2C — 0.2411Ji and J2L — 0.5206Ji respectively. Though no real 
phase transition occurs for J2 > J2C in the model, the dimerized phase in 
the region J2 > J2L can be distinguished from that at J2C < J2 < J2L 
due to the existence of this incommensurability. In this respect, the sponta- 
neously dimerized phase for J2 > J2L with deconfined massive spinous and 
incommensurate correlations represents a remarkable ID spin liquid phase 
stabilized by frustration. 

This incommensurability induced by frustration is, in fact, quite general 
in ID and not restricted to the spin-1/2 case. In particular, this phenomenon 
appears also in the spin-1 case. The zero-temperature phase diagram of the 
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spin-1 Ji — J2 Hciscnberg chain (|24|) has been investigated numerically 
0-92 158 193 194 195 1961 ^^-^^ ^ spectral gap exists fo r all finite value of J2. 
A first-order transition occurs at i'h.l J ])t .^.Q ■ 7M separating a Hal- 
dane phase from a double Haldane phase^ ^ "^ 1 ^ " v l 1 y 9l xhe quantum number 
carried by elementary excitations is 5* = 1 in full agreement with the semi- 
classical approach described above. Frustration plays its trick here by 
giving rise to a similar onset of incommensurability as in the spin-1/2 case. 
A disorder point of the first kind occ urs a t {J2I Ji)d — 0.284 together with 
a Lifshitz point (J2/Ji)l — 0.3725 after which the static structure 
factor develops a two-peak structure. 

3.3. Field theory of spin liquid with incommensurate 
correlations 

Frustration represents thus a novel mechanism in one dimension for gen- 
erating incommens urab ility as external magnetic fields or Dzyaloshinskii- 
Moriya interaction xhe zigzag ladder representation (Fig. of the 
Ji — J2 model enables one to investigate, by a weak coupling approach 
Ji <C J2, the main characteristics of frustration and to shed light on the 
mechanism that gives rise to incommensurate correlations in the large J2 
regime. In the spin-1/2 case, the onset of incommensurability and the pres- 
ence of deconfined massive spinons might be understood starting from the 
limit where the two spin-1/2 AF Heisenberg chains are decoupled. 

The continuum lim it of thejS = 1/2 two-leg ladder has been anal- 

ysed by several ^.on^M^^^^^^, 

The interacting part of 
the Hamiltonian density of the low-energy field theory reads as follows: 

Hint ^ 91 i^JiL ■ J2R + J2L ■ JiRj +g2ni- dxn2, (38) 

where Ja.L.R (a = 1,2) are the left-right su(2)i currents corresponding to 
the continuum description of the ath spin-1/2 AF Heisenberg chain. In Eq. 
(I38II ■ Ha denotes the staggered magnetization of the spin density of the 
chain with index a — 1,2. In the continuum limit, this fia field identifies 
with the vector part of the primary field, with scaling dimension A„ = 1/2, 
of the su(2)i WZNWCFT transforming according to the fundamental rep- 
resentation of SU(2)^EnD4 10_ ^Yiis continuum description, frustration 
suppresses geometrically the standard backscattering contribution Hi ■ H2 
which governs the low-energy physics of the two-leg spin ladder. This term 
represents a strongly relevant perturbation of scaling dimension 1. It corre- 
sponds to the energy operators of the underlying Ising models of the low- 
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energy description (|16|l of the two-leg spin ladderlSHE21 As seen in Section 
2, this backscattering contribution leads to the formation of a spin gap in 
this model as well as the confinement of the original massless spinons into 
optical magnons. In presence of frustration, the interacting part (|38|l is now 
only marginal relevant and consists of two terms of different nature. The 
first one, with coupling constant gi, is a current-current interaction similar 
to that which appears in the effective field theory H29|l of the spin- 1/2 J1 — J2 
Heisenberg chain in the weak coupling lim it Jo <C Ji ■ The second marginal 
contribution in Eq. H38|l . called twist terml ^^ l, Ctwist = ni ■ dxn2 is a novel 
parity-breaking perturbation which contains the staggered magnetizations 
of each chain and also a spatial derivative. In contrast to the current-current 
perturbation, this twist term is not a scalar under the Lorentz transforma- 
tion in 1 -f 1 dimensions but behaves as a vector. In the CFT jargon, this 
kind of perturbation is characterized by a non-zero conformal spin 5 = ±1 
together with its scaling dimension A = 2. In fact, this twist perturbation 
has been forgotten in the first bosonization analysis of the two-leg zigzag 
spin ladder ^'^^ ^"iQ It ha s been discovered by Nersesyan et al. ^'^'^ and in- 
dependently by AUenl^SI, Jt is interesting to observe that a similar operator 
appears in the effective field theory approach to the spin-1 two-leg zigzag 
ladder in the large J2 limit ^ . This twist perturbation is thus the hall- 
mark of frustration in the low-energy description of frustrated spin ladders 
with zigzag interchain interaction. 

The effect of such a non-zero conformal spin perturbation is non-trivial 
since the usual irrelevant versus relevant criterion of perturbative field the- 
ory does not hold for s uch a non-scalar contribution (see for instance the 
discussion in the book 1^1). In the case of a 5* = ±1 perturbation, it is 
expected that the generic effect of this contribution is the onset of in- 
commensurability. In this respect, a simple example of incommensurability 
arising from such a perturbation is the spin- 1/2 XXZ Heisenberg chain in 
a magnetic field along the z-axis. The Hamiltonian density of th is model, 
obtained within the bosonization approach, reads as follows 



V 

^-2 



{d.,<^f + {dx<df -5Cos/3$-/i5x$. (39) 



Here the S" = ±1 perturbation is described by the uniform part of the 
spin density dx'^. It is well known that in this model, when the magnetic 
field h is increased, a commensurate-incommensurate phase transition takes 
place with the appearance of an incommensurate phase with critical cor- 
relation functions 205 206 207 208 x}jjg transition occurs for a finite or 
vanishing magnetic field depending on the relevance or irrelevance of the 
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cos/3$ term. A similar incommensurability generated by a S" = ±1 confor- 
mal spin perturbation o ccurs in the spin-1 /2 XXZ Heisenberg chain in a 
transverse magnetic field l^^^^^^^^ffl jn g, model of two spinl ess Luttinger 
chains weakly-coupled by a single-particle interchain hopping l^^E^^EEl 
and in the quantum axial next-to-nearest neighbor Ising (ANNNI) chain in 
a transverse magnetic field l^i^I Finally, it is important to note that there 
exists some exactly solvable CFT models perturbed by = ±1 confor- 
mal spin term for which the presence of incommensurability can be shown 
non-perturbatively "^^"^ ^Ib 

All these results suggest tha t the twist term Otwist ^ ni ■ dxn2, as 
proposed by Nersesyan et al. should be at the origin of the incommen- 
surability found numerically in the spin- 1/2 Ji — J2 Heisenberg chain (|24l) 
when J2 > O.52J1. To this end, this twist term and the current-current 
interaction of Eq. H38I) can be expressed in terms of the four Majorana 
fermions Cal'*^ ~ 0,1,2,3 of the continuum limit l|16|l of the standard 
two-leg spin ladder | 1^8 | 199 | jj^ particular, the current-current interaction 
is built from of the following tensor: O^^ = ^r^Ur^I with a, 5 = 0, 1, 2, 3 
and a ^ b. The twist perturbation is also local in terms of these Majorana 
fermions but with a different structure. A typical term that enters its ex- 
pression is O^tf^t = ^r^l^UI with a, 6, c, d = 0, 1, 2, 3 and a ^ 6 ^ c ^ d. 
The non-zero conformal spin of this twist term is reflected here by the dif- 
ferent number of right and left fermions in this expression. The field theory 
(I38|l . expressed in terms of these four Majorana fermions, turns out not to 
be integrable. The renormalization group (RG) flow analysis reveals that 
the current-current interaction and the twist term are equally important in 

the IR limit. They reach a strong coupling regime simultaneously with a 
1 Q8 901 909 

fixed ratio ^"^^ ^'J^. The nature of this strong coupling regime is still an 
open problem since the emerging IR field theory is not integrable and one 
cannot disentangle the effect of the two interactions of Eq. (|38|l . However, 
it is very tempting to explain the main characteristics of the spontaneously 
dimerized phase with incommensurate correlation of the spin-1/2 Ji — J2 
Heisenberg chain in the large J2 limit from these two contributions. On 
one hand, the current-current contribution of Eq. (jSHJ, equivalent to two 
decoupled SU(2) Thirring model, has similar properties as the field theory 
(|29ll obtained in the weak coupling regime J2 <C Ji- The spontaneously 
dimerization and the existence of massive deconfined spinons should re- 
sult from this interaction. In fact, this current-current interaction appears 
alone, without a twist term, in the continuum limit of a frustr ated t wo-leg 
spin ladder with crossings along a special line of the couplings 1^^. It has 
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199 



been shown that this model, in the weak coupling limit, is characterized 
by a weak spon taneously dimerization with massive spinous as elementary 
excitations'^"^. Finally, as already stressed, the twist perturbation of Eq. 
should be at the origin of the incommensurability in the large J2 limit 
In this respect, the effect of the twist term of the two- leg zigzag lad- 
der has been analysed within a RPA approach and it leads indee d to some 
incommesurability behavior in the spin-spin correlation function l "^^^! 

A simple way to disentangle the effects of the current-current interaction 
and the twist perturbation is to introduce an exchange anis qtrop y which 
makes the twist contribution more relevant in the RG sense l -^^^^ l The ex- 
treme case is the XY version of the spin- 1/2 Ji — J 2 Heisenberg chain with 
Hamiltonian: 



n n 

In the ladder limit when J2 S> Ji, the low-energy physics of this model 
can be analysed by means of the bosonization approach. Introducing two 
bosonic fields th e resulting bosonic Hamiltonian density of the model 
reads as follows' I 

HxYC^lY. + idxOaf) + g d,e+ sin (\/2^e-) , (41) 

where Q± are the dual fields associated to the bosonic fields The Hamil- 
tonian (|41l) describes a nontrivial field theory with a relevant 5 = ±1 con- 
formal spin twist perturbation with scaling dimension Ag = 3/2. In this 
strong anisotropic XY case, the current-current perturbation is less domi- 
nant and can be safely neglected to derive the nature of the phase of the 
model when J2 3> Ji- The presence of incommensurability in the system 
can then be found using a mean-field analysis of the model by decoupling 
the two pieces of the twist term of Eq. (|41() . In particular, the leading 
asymptotic behavior of the transverse spin-spin correlation functions of the 
model obtained by this mean-field approach is given 

byim 

(5t(x)5-(0))^^^|^, a =1,2, (42) 

which displays an incommensurate critical behavior with an oscillating fac- 
tor: q — n/ao ^ (Ji/J2)^. The physical picture that emerges from this 
mean-field analysis is the existence of a critical spin nematic phase that 
preserves the U(l) and time-reversal symmetries and displays long-range 
chiral ordering in its ground state: ((^n+i A Sn)z) 7^ 0. It is important to 
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note that this spin nematic phase does not break the time-reversal sym- 
metry but spontaneously breaks a Z2 symmetry of the model which is a 
tensor product of a site-parity and link-parity symmetries on the two chains: 
P^^^ X -P^^"*. The z-component of the spin current J^^ associated to the ath 
spin-1/2 XY chain (a = 1,2) takes a non -zero expectation value in the 
ground state of the mean field Hamiltonian 11^^ 

{JD = (JL) = {d.Q+) + 0. (43) 

As a result, this produces a picture of local nonzero spin currents polar- 
ized along the z-anistropy axis circulating around the triangular plaquettes 
of the two-leg zigza g spin ladder. This phase has been found numerically 
by Nishiyama I ^ 1 who has investigated the existence of chiral order of 
the Josephson ladd er with half a flux quantum per plaquette, and also by 
Hikihara et al. I^EI by means of the DMRG approach. This critical incom- 
mensurate pha se is the quantum analogue of the classical spiral phase with 
chiral ordering It is thus natural to expect that this phase is not re- 
stricted to the spin-1/2 case but should exist in the range of parameters of 
the model (|^ in the general spin case. Indeed, the existence of this criti- 
cal spin nem atic phas e in the spin-S' case has been show n num erically for 
S" = 1, 3/2, 2EIlE2ni^ by means of a semiclassical method and finally 
by the Abelian bosonization approach ' ' of a general spin S introduced 
by SchulzEi. 

3.4. Extended criticality stabilized by frustration 

In addition to incommensurability, a second characteristic of classical 
canted magnets is the presence of three gapless spin- wave modes instead of 
two as in the Neel colinear state. An interesting question is whether frustra- 
tion can lead to new type of emerging quantum criticality not encountered 
in unfrustrated spin chains or ladders. In this last case, the critical behav- 
ior of half-integer AF Heisenberg spin chains or odd-legged spin ladders 
is characterized by one gapless bosonic mode i. e. by a CFT with central 
charge c = 1. As it has already been pointed out in the previous section, 
a striking effect of frustration in the continuum limit stems from the fact 
that the low-energy effective field theory is mainly governed by marginal 
interactions. New type of IR critical behaviors may result from the delicate 
balance between these marginal perturbations. In the following, we shall 
give some examples of critical phases and quantum critical points with ex- 
tended criticality stabilized by frustration in ID spin systems. 
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3.4.1. Critical phases with SU(N) quantum criticality 

A first example of a system with extended criticality is the ID spin-orbital 
model (I21|l in a regime where frustration reveals itself K J. At the special 
point J = K/A, the Hamiltonian H21|l can be expressed in terms of a product 
of two-body permutation operator in Si and S2 subspaces: 

= jT.pSxt"'pSxt"'- (44) 

i 

Since this Hamiltonian exchanges both 5*1..; and 6*2. i spins at the same time, 
the spin-orbital model at J = K/4: is not only SU(2) x SU(2) symmetric 
but actually has an elarged SU(4) symmetry which unifies the spin and 
orbital degrees of freedom l^^vl^^'^ l ^^-.'l More precisely, the Hamiltonian 
(|44|l can be recasted as an AF Heisenberg spin chain with SU(4) spins up 
to a constant: 

15 

Hsu(4) ^^EE^'^^+i' (45) 

i A=l 

where are the 15 generators belonging to the fundamental represent ation 
of SU(4). This model is exactly solvable by means of the Bethe ansatzE^ 
and its low-energy spectrum consists of three gapless spinous with wave 
vectors ±7r/2ao and tt/oo '^'^'^ . As shown by Affleck 2^^, the critical 
theory corresponds to the su(4)i WZNW model with central charge c = 
3 (three massless bosonic modes). The existence of this quantum critical 
point with a SU(4) symmetry allows us to study the spin-orbital model 
(|21|l by a continuum description in an intermediate coupling regime J ^ 
K where frustration shows off. In this respect, it is im portant to notice 
that, according to the Zamolodchikov's c theorem |229l ^ the SU(4) critical 
point with central charge c = 3 cannot be reached by a RG trajectory 
starting from the decoupling limit {K ~ 0) of two spin-1/2 AF Heisenberg 
chains with total central charge c = 2. Stated differently, the physics in 
the neighborhood of the SU(4) point cannot be understood in terms of 
weakly coupled S* = 1/2 Heisenberg chains. The strategy to tackle with this 
intermediate coupling regime is to start from the SU(4) Hubbard chain at 
quater-fiUing and apply the bosonization approach to obtain the continuum 
description of the spin densities 6^1,2(2;) at the SU(4) point in the large 
Coulomb repulsion limit. The low-energy field theory which describes small 
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deviations from the SU(4) symmetric point can then be derived [ ^^^ l ^'^l l 
The resulting effective Hamiltonian density associated to the symmetry 
breaking scheme SU(4) SU(2) x SU(2) of the spin-orbital model ^ 
reads as follows E^Hl 



'HeS = ^ (^Cr ■ dxisR. — isL ' dx^shj ^ {^oR ' d^^oR — ^oL ' dx^oL 

^^isR ■ ishj + 92 (Cofi • Col) + 53 (Csfl ■ ClJ (^ofl ■ ?ol) , (46) 



+ 51 

where ^sR,l and ^oR,l are two triplet of Majorana fermions which act re- 
spectively in the spin and orbital sectors. In Eq. (|46|l . we have considered a 
more general situation than the model H21|) by allowing the exchange in spin 
( Ji) and orbital (J2) channels to be different: the three coupling constants 
gi are independent. The interaction of the low-energy field theory H46|l is 
marginal and describes two SO (3) Gross-Neveu^^ models marginally cou- 
pled. This field theory is not integrable so that one has to recourse to per- 
turbation theory to elucidate the phase diagram of the spin-orbital model 
in the vicinity of the SU(4) symmetric point. In fact, the all-order beta 
functions of the fi eld th eory H46|l can be determined using the approach 
of Gerganov et al. ES21 These authors have computed the all-order beta 
functions of a general model with anisotropic current-current interactions 
in a special minimal scheme prescription. The model H46|) can be viewed 
as a current-current interaction corresponding to the symmetry breaking 
scheme S0(6 ) — » S0(3)x SO (3). The application of the general formula 
given in Ref. ^ confirms the conclusions derived from the one-loop cal- 
culation '^'^^ '^'^^ which reveals the existence of two phases with different 
remarkable properties. 

A first phase, for Ji ~ J2 > K/i, has a spectral gap and the ground 
state has a similar staggered dimerization as the weak-coupling phase of Fig. 
0] Provided the anisotropy is not too large Ji 7^ J2 in the vicinity of the 
SU(4) symmetric point, the massive phase displays an approximate S0(6) 
~ SU(4) enlarged symmetry. Indeed, by neglecting the velocity anisotropy 
in Eq. (|46|l . the RG equations reveal a fiow to strong coupling with an 
attraction along a special directio n wi th enlarged symmetry described by 
the following Hamiltonian density ESDI 

6 / ^ \ ^ 

HiR^-^E(^Sa.S,-e£a,e£) + 5JE^«^M ' (47) 

a=l \a=l / 

with 5* > 0. This Hamiltonian is identified as the SO (6) Gross-Neveu model 
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which is a massive integrable field theory. Its spectrum is known and con- 
sists of the fundamental fermion, with mass A/, together with a kink and 
anti-kink with mass mkink = M/\/2 . The initial model (gHll, SU(2) x 
SU(2) symmetric, acquieres thus in the IR limit an enlarged SO (6) sym- 
metry. In more physical terms, the spin and orbital degrees of freedom are 
unified and described by a same multiplet with six components. A similar 
example of symmetry restauration by interactions is the emergence of a 
SO (8) symmetry in weakly-coupled two-leg Hubbard ladder at half-filling 
1233 g^^^ tj^g su(4) Hubbard chain at half-filfing 231 

The second phase with Ji ~ J2 < K/A has striking properties. All cou- 
plings in Eq. I|46|l flow to zero in the IR limit and the interaction is marginal 
irrelevant. The six Majorana fermions are thus massless and the phase dis- 
plays extended quantum criticality characterized by a central charge c — Z. 
Spin-spin correlation functions decay algebraically with exponent 3/2 and 
exhibit a four-site periodicity {2kp — 7r/2ao). The critical behavior at the 

1.5 
1 

0.5 







-0.5 



-1 -0.5 0.5 1 1.5 

X 

Fig. 11. Phase diagram of the spin-orbital model obtained by DMRG; x = J2/K and 
y = Ji/K; taken from Ref. 231. 

SU(4) point extends to a finite region of the phase diagram of the spin- 
orbital model. In this respect, this remarkable gapless phase with extended 
quantum criticality c — 3 represents a new universality class, stabilized by 
frustration, in spin chains and spin ladders. It is worth noting that this 
phase has a maximum of gapless modes allowed by the classical structure 
of a spiral (three Goldstone modes). 

The phase diagram of the anisotropic spin-orbital model H21|l with 
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Ji ^ J2 has been investigated numerically by means of the DMRG tech- 
nique ^^'^ ^"^^ '^'^-'_jFiKmeE] represents the resulting T — phase diagram 
obtained in Refs. I '^^ 1 I In phase I, both spin and orbital degrees of free- 
dom are in fully polarized ferromagnetic states. In phase II, the orbital 
degres of freedom are still in a ferromagnetic state whereas the spin de- 
grees of freedom are now critical (criticality of a spin- 1/2 AF Heisenberg 
chain) and vice versa in phase III. Phase IV corresponds to the staggered 
dimerized phase of Fig. 0| while phase V is the c = 3 gapless phase with 
extended quantum criticality. In fact, this critical phase has been first dis- 
covered only along a special line —K/A < Ji = J2 < in the DMRG 
calculation of Pati et al. 

11321. Additional DMRG worksESHSIl 

have shown 

the existence of this gapless phase in an extended region (Phase V) of the 
phase diagram (Fig. Illfl. 

A second model which displays extended quantum criticality induced by 
frustration is the spin-1 bilinear-biquadratic chain @ with a biquadratic 
coupling constant f3 > 1. In this regime, frustration manifests itself and the 
classical ground state is an incommensurate spiral. The Haldane phase of 
the spin-1 bilinear-biquadratic chain @ with \P\ < 1 e nds at th e integrable 
point with /? = 1, the so-called Uimin-Lai-SutherlandESSESEI point. This 
model can be expressed in terms of SU(3) spins and takes the form of a 
SU(3) AF Heisenberg spin chain: 

Hsu(3)- JEE^^^^I' (48) 

i A=l 

where are the 8 generators belonging to the fundamental representa- 
tion of S U(3). This model is exactly solvable by means of the Bethe ansatz 
l ^^?(j|^^t>l and its low-energy spectrum consists of two gapless spinous with 
wave vectors ±27r/3ao As shown by Affleck ^^^^ the critical theory 
corresponds to the su(3)i WZNW model with central charge c = 2 (two 
massless bosonic modes) . The existence of this quantum critical point with 
a SU(3) symmetry allows to study the spin-1 bilinear-biquadratic chai n Q 
in the vicinity of /3 ~ 1 by means of a field theory approach ESZI xhe 
low-energy effective field theory is described by marginal current-current 
inter actions associated to the symmetry breaking scheme SU(3) SU(2) 
12^. The one-loop RG flow near /3 ~ 1 reveals the existence of two distinct 
phases. On one hand, a phase with a dynamical mass generation when 
(3 < 1 which signals the onset of the Haldane phase. On the other hand, 
a massless phase for /3 > 1 where interactions are marginal irrelevant. The 
SU(3) quantum criticality of the Uimin-Lai-Sutherland point, with two gap- 
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less bosonic modes, extends to a finite region of the phase diagram of the 
generalized spin-1 Heisenberg chain 0. This massless phase has been first 
pred icted in a numerical investigation of this model by Fath and Solyom 
I I This massless phase with approximate SU(3) symmetry, stabilized by 
frustration, is similar in spirit to the previous SU(4) massless phase of the 
spin-orbital model 121|) . Finally, it is worth noting that this c = 2 massless 
phase appears also in the phase diagram of the so-called spin tube model, 
a three-leg spin ladder with frustrated periodic interchain interaction, in a 
magnetic field 



3.4.2. Chirally stabilized critical spin liquid 

Frustration can also induce an exotic quantum critical point in frustrated 
spin ladders. A possible lattice realization of this phenomenon is a model 
(see Fig. IT^ of three = 1/2 AF Heisenberg spin chains weakl y cou pled 
by on- rung J± and plaquette-diagonal J± interchain interactions'^^. The 
Hamiltonian of this model reads as follows: 



^11 Sa^i • + Jl_ 52,1 • i^SiA + Sz^i 



i a—1 



Jx J2 [{^^'- + ^^ i) ■ ^2,^+1 + {Si,,+ i + 53,.+ l) • S2,^] . (49) 



In the continuum limit with J±, Jx <C Jii, the Hamiltonian of the lattice 




Fig. 12. Three-leg spin ladder with crossings. 
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model H49f) takes the form: 



Zttv 



ZTTV 's—^ / -h -h \ ^ 

— l^-^afl + JaL) +9^2- («! + na) 

a— 1 



(50) 

with g = ao{J± + 2Jx) and g — uq (J± — 2Jx)- In Eq. jSUJ, JaR,L 
(a = 1,2,3) are the right and left chiral su(2)i currents which accounts 
for the low-energy description of the uniform part of the spin density of 
the ath spin- 1/2 chain. It is interesting to note that in Eq. (|50|l . there is no 
marginally relevant twist perturbation nidxn2 which appears in the contin- 
uum description of the two-leg spin ladder with a small zigzag interchain 
coupling (|38|l . The two interaction terms in Eq. (|50|l are of different nature. 
On one hand, the first contribution, with coupling constant g, is a strongly 
relevant perturbation with scaling dimension Ag = 1. On the other hand, 
the second term is a current-current interaction which is only marginal and, 
as long as g is not too small, can be discarded. As a result, for generic val- 
ues of g and the low-energy physics of the model H49(l will be essentially 
that of the standard three-leg ladder, and frustration will play no role (ex- 
cept for renormalization of mass gaps and velocities) . The important point 
here is that in contrast with non-frustrated ladders, the two coupling con- 
stants g^g can vary independently, and there exists a vicinity of the line 
J_L = 2 Jx {g = 0) where the low-energy properties of the model are mainly 
determined by current-current interchain interaction. 

Along the special lin e J i = 2 Jx , the low-energy effective Hamiltonian 
(150(1 simplifies as follows 



j J I 

Vim 



g J2R ■ II + J2L ■ Ir 



2L ; 



(51) 



+ T 

where ^ are Majorana fermions which are associated to the Z2 (1 — > 3) 
discrete interchange symmetry between the surface chains labelled a = 1 , 3 
in the lattice model H49|) : the total chiral current of the surface chains is 
noted Ir^l = Jir,l + J^rx and corresponds to a su(2)2 WZNW current. 
The Hamiltonian (|51|l has an interesting structure. First, the Majorana 
fermions ^ do not participate in the interaction and remain thus critical 
in the IR limit. These massless degrees of freedom can be interpreted as an 
effective two-dimensional Ising model at T = Tc which accounts for singlet 
excitations between the surface chains. All non-trivial physics of the model 
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is incorporated in the current-current dependent part of the Hamiltonian 
(|51|l . denoted by TCcc, which describes marginaUy coupled su(2)i and su(2)2 
WZNW models. Moreover, this Hamiltonian separates into two commuting 
and chirally asymmetric parts: Ticc = T^i +'^2, (['^i,'W2] = 0), where 

-Hi^ — Pn + — -^L+g Ir ■ J2L, (52) 

and Ti.2 is obtained from Tii by inverting chiralities of all the curr ents. The 
model (152(1 is integrable by means of the Bethe-ansatz approach ^^^^21 
A simple RG analysis shows that, at g > 0, the interaction is marginally 
relevant. Usually the development of a strong coupling regime is accompa- 
nied by a dynamical mass generation and the loss of conformal invariance 
in the strong coupling limit. However, here due to the chiral asymmetry of 
Til , it turns out that the effective interaction flows towards an intermediate 
fixed point where conformal invariance is recovered with a smaller central 
charge. This critical behavior has been identified as the universality class 
of chirally stabilized fluids, introduced by Andrei, Douglas, and Jerez 1^^. 
In the case of the model H52|) . the symmetry of the IR fixed point, obtained 
from the Bethe-ansatz analysis B*^ , turns out to be su(2)i |r x Z2\l- This 
result can also be derived using a Toulouse point approach* As a whole, 
taking into account of the contribution of the Majorana fermions ^, the 
model (|51ll displays critical properties characterized by a fixed point with 
a su(2)i X Z2 X Z2 symmetry and a c = 2 central charge. This emerging 
quantum criticality can be interpreted as the criticality resulting from an 
effective 5 = 1/2 AF Heisenberg chain and two decoupled critical Ising 
models. Frustration, introduced by the diagonal exchange interaction Jx , 
shows up in a non-trivial way by stabilizing a quantum critical point with 
central charge c — 2. This IR behavior differs from the standard c = 1 
quantum criticality in unfrustrated three-leg spin ladder by the two gapless 
non-magnetic, singlet, degrees of freedom described by the Ising models. 
The physical properties of the model at this c = 2 critical point can be de- 
termined by a Toulouse point approach of the model 1(5 l|l The slowest 
spin-spin correlation functions of the model at the IR critical point corre- 
spond to the staggered correlations between the spins of the surface chains 
which decay with a universal exponent 3/2. As a consequence, the low- 
temperature dependence of the NMR relaxation rate scales as 1 /Tj ~ -s/T 
in contrast to 1/Ti ^ const for the spin-1/2 AF Heisenberg chain The 
c = 2 quantum critical point, induced by frustration, describes thus a new 
universality class in spin ladders. In addition, the Toulouse point analysis 
enables us to investigate the effect of the neglected backscattering term 
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^2 ■ (^1 + ^3) of the model (|5U|I in the vicinity of the hne g = 0. The main 
effect of this operator is to open a spectral gap in the Ising degrees of free- 
dom of the c = 2 fixed point but has no dramatic effect on the effective 
5 = 1/2 AF Heisenberg chain. The c = 2 IR fixed point is thus unsta- 
ble with respect to the interchain backscattering term and the model (150(1 
will display IR critical properties governed by the c = 1 fixed point of the 
three- leg spin ladder or the 5 = 1/2 AF Heisenberg spin chain. 

Finally, it is worth noting that the zero-temperature phase diagram of 
the three-leg spin ladder with cro ssing s has been recently investigated 
by means of the DMRG approachl I Two critical phases have been found 
corresponding to the 5 = 3/2 AF Heisenberg chain and the three-leg spin 
ladder. Frustration, introduced through the diagonal interchain interaction 
Jx , induces a quantum phase transition between these two critical phases 
belonging to the standard c = 1 universality class. At large Jx , the transi- 
tion has been found to be of first order. In the weak coupling regime, the 
transition occurs at J± = 2Jx when Jx < 0.8J|| in full agreement with the 
field theoretical description (|50|l of the model. Extended DMRG calcula- 
tions are required to fully characterize the nature of the phase boundary 
at J_L = 2 Jx in the weak coupling regime in particular to verify the chiral 
spin liquid behavior with c = 2 criticality proposed in Ref. -2411 

4. Concluding remarks 

The suppression of magnetism by quantum fluctuations in one dimension 
gives rise to a large variety of zero-temperature spin liquid phases with 
striking different physical properties. The canonical examples are the emerg- 
ing quantum criticality of the spin-1/2 AF Heisenberg chain, controlled by 
gapless elementary excitations with fractional quantum numbers (spinous), 
and the formation of an incompressible Haldane spin liquid phase in the 
spin-1 case with optical 5=1 magnon excitations. Several other types of 
incompressible spin liquid phases can be stabilized with different dynami- 
cal or optical properties. A distinct spin liquid behavior may be revealed 
in the form of the dynamical structure factor with the presence or not 
of a sharp spectral peak which directly probes the elementary nature of 
the triplet excitation of the phase. In this respect, the staggered dimerized 
phase of the two-leg spin ladder with a biquadratic interchain exchange 
represents a non-Haldane spin liquid phase even though a spectral gap is 
formed by quantum fluctuations. The main difference between this phase 
and the Haldane phase of the spin-1 Heisenberg chain stems from the com- 
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posite nature of the spin-flip S — I excitation due to the presence of a 
two-particle threshold in the dynamical structure factor instead of a sharp 
magnon peak. Two gapped spin liquid phases can also be distinguished at 
the level of the topology of the short-range valence bond description of their 
ground states. A spin liquid phase may exhibit a ground-state degeneracy 
depending on the nature of boundary conditions used. In particular, for 
open boundary conditions, chain-end spin excitations can result from this 
ground-state degeneracy leading to well defined satellite peaks in the NMR 
profile of the system doped with non-magnetic impurities. An example of 
this topological distinction between two gapped phases is provided by the 
two-leg spin ladder. In the Haldane phase of the open two-leg ladder with 
a ferromagnetic rung interchain Jl < 0, = 1/2 edge states are formed 
whereas these chain-end excitations disappear in the rung singlet phase of 
the ladder with J± > 0. 

Frustration, i.e. the impossibility to satisfy simultaneously every pair- 
wise interaction, can induce additional types of spin liquid phases with 
exotic properties. In the spin-1/2 case, one of the most striking effect of 
frustration is the stabilization of a spontaneously dimerized phase with 
massive deconfined spinon excitations. Frustration represents thus a direct 
route to fractionalization. The existence of these deconfined spinon excita- 
tions has important consequences, as for instance, in the study of doping 
effects of such a spin liquid phase with non-magnetic impurities. On general 
grounds, it is expected that no free spin degrees of freedom are generated 
arou nd non-magnetic impurities due to the presence of deconfined spinous 
In the case of the spontaneously dimerized phase of the spin-1/2 Ji- 
J2 Heisenberg chain, the absence of induced free spin degrees of freedom 
or edge states has been shown recently A second remarkable effect of 
frustration is the onset of incommensurability in spin chains whose classical 
ground state has a spiral structure. In this respect, frustration is a novel 
mechanism in one dimension for generating incommensurate behavior as 
external magnetic fields or Dzyaloshinskii-Moriya interaction. An example 
of such a spin liquid phase with incommensurate correlation is the quantum 
Ji — J2 Heisenberg chain in the large next-nearest neighbor limit. When 
the SU(2) symmetry of this model is explicitely broken, an incommensurate 
phase is produced by frustration with local non-zero spin currents, polar- 
ized along the anisotropic axis, circulating around triangular plaquettes. 
This phase is the quantum signature of the classical spiral phase induced 
by frustration. Finally, a last effect of frustration, discussed in this review, 
is the possible realization of a new type of emerging quantum criticality in 
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spin chains and ladders. In the continuum description, the low-energy effec- 
tive field theory of a frustrated spin chain is mainly governed by marginal 
interactions. An extended IR critical behavior characterized by a CFT with 
central charge c > 1 may result from the delicate balance between these 
marginal contributions. 

Regarding perspectives, a natural question to raise is the existence of 
a novel ID spin liquid phase with unbroken SU(2) spin symmetry which 
displays physical properties not described in this review. In particular, an 
interesting possibility is the realization of a ID version of the chiral spin 
liquid phase ESI which breaks spontaneously the time-reversal symmetry. 
In fact, such a phase has been identified recently in the phase diagram of 
the two-leg spin ladder with a four-spin cyclic exchange »246t247, jj^ ^j^g 
large ring-exchange limit, a spin liquid phase characterized by a non-zero 
scalar-chirality operator {Si^i ■ (5*2, i A Si^i+i)), breaking both p arity and 
time-reversal symmetries, has been found in DMRG calculations l^46 | 2 47| 
Using duality arguments, the authors of Ref. have also exhibited sev- 
eral spin ladder models with exact ground-state which display spontaneous 
breaking of the time-reversal symmetry. Finally, a central issue is the possi- 
bility to stabilize a two-dimensional SU(2)-invariant spin liquid phase with 
exotic properties starting from the one-dimensional limit. Recently, a model 
of two-dimensional frustrated spin system with SU(2) spin symmetry and 
strong spatially anisotropy has been introduced by Nersesyan and Tsvelik 
This model can also be interpreted as a collection of weakly coupled 
spin chains. It has been shown that, in the limit of infinite number of 
chains, this system displays a spin liqui d ph ase with massive deconfined 
fractional spin-1/2 excitations (spinous) l ^^*^ ! The existence of spinous in 
two- dimensional spin systems has also been found in a crossed-chain model 
which is a 2D anisotropic version of the p Yrochlore lattice. Remark- 
ably enough, for a range of coupling constants l ^49 | 25U I^ ^^lis system is a 
non-dimerized spin liquid with deconfined gapless spinous as elementary 
excitations. This spin liquid phase, stabilized by frustratio n, is an exam- 
ple of a SU(2) version of a sliding Luttingcr liquid phase I I We hope 
that other two-dimensional spin liquid phases with exotic properties will 
be reported in the near future. 
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